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METHODS AND COMPOSITIONS FOR DETECTION OR 
QUANTIFICATION OF NUCLEIC ACID SPECIES 



HROSS REFERENCE TO RELATED APPLICATIONS. 

5 This patent application is a continuation-in-part of 

U.S. Patent application Serial No. 08/812,951, filed on March 
4, 1997, and a continuation-in-part of U.S. Patent application 
Serial No. 08/784,747, filed on January 16, 1997. 

10 FIELD OF THE INVENTION 

This invention relates in general to methods and 
apparatus for nucleic acid analysis, and, in particular, to 
methods and apparati for nucleic acid analysis. 

15 BACKGROUND 

The rate of determining the sequence of the four 
nucleotides in nucleic acid samples is a major technical 
obstacle for further advancement of molecular biology, 
medicine, and biotechnology. Nucleic acid sequencing methods 
20 which involve separation of nucleic acid molecules in a gel 
have been in use since 1978 . The other proven method for 
sequencing nucleic acids is sequencing by hybridization (SBH) . 

The traditional method of determining a sequence of 
nucleotides { i.e., the order of the A, G, C and T 
25 nucleotides in a sample) is performed by preparing a mixture 
of randomly- terminated, differentially labelled nucleic acid 
fragments by degradation at specific nucleotides, or by 
dideoxy chain termination of replicating strands. Resulting 
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nucleic acid fragments in the range of 1 to 500 bp are then 
separated on a gel to produce a ladder of bands wherein the 
adjacent samples differ in length by one nucleotide. 

The array-based approach of SBH does not require 
5 single base resolution in separation, degradation, synthesis 
or imaging of a nucleic acid molecule. Using mismatch 
discriminative hybridization of short oligonucleotides K bases 
in length, lists of constituent K-mer oligonucleotides may be 
determined for target nucleic acid- Sequence for the target 
10 nucleic acid may be assembled by uniquely overlapping scored 
oligonucleotides . 

There are several approaches available to achieve 
sequencing by hybridization* In a process called SBH Format 1, 
nucleic acid samples are arrayed, and labeled probes are 

15 hybridized with the samples. Replica membranes with the same 
sets of sample nucleic acids may be used for parallel scoring 
of several probes and/or probes may be multiplexed. Nucleic 
acid samples may be arrayed and hybridized on nylon membranes 
or other suitable supports . Each membrane array may be reused 

20 many times. Format 1 is especially efficient for batch 
processing large numbers of samples . 

In SBH Format 2, probes are arrayed at locations on 
a subsrate which correspond to their respective sequences, and 
a labelled nucleic acid sample fragment is hybridized to the 
25 arrayed probes. In this case, sequence information about a 
fragment may be determined in a simultaneous hybridization 
reaction with all of the arrayed probes. For sequencing other 
nucleic acid fragments, the same oligonucleotide array may be 
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reused. The arrays may be produced by spotting or by in situ 
synthesis of probes. 

In Format 3 SBH, two sets of probes are used. In 
one embodiment , a set may be in the form of arrays of probes 
5 with known positions, and another, labelled set may be stored 
in multiwell plates. In this case, target nucleic acid need 
not be labelled. Target nucleic acid and one or more labelled 
probes are added to the arrayed sets of probes. If one 
attached probe and one labelled probe both hybridize 

10 contiguously on the target nucleic acid, they are covalently 

ligated, producing a detected sequence equal to the sum of the 
length of the ligated probes. The process allows for 
sequencing long nucleic acid fragments, e.g. a complete 
bacterial genome, without nucleic acid subcloning in smaller 

15 pieces. 

In the present invention, SBH is applied to the 
efficient identification and sequencing of one or more nucleic 
acid samples. The procedure has many applications in nucleic 
acid diagnostics, forensics, and gene mapping. It also may be 
20 used to identify mutations responsible for genetic disorders 
and other traits, to assess biodiversity and to produce many 
other types of data dependent on nucleic acid sequence. 

SUMMARY OF THE INVENTION 

25 The present invention provides a method for 

detecting a target nucleic acid species including the steps of 
providing an array of probes affixed to a substrate and a 
plurality of labeled probes wherein each labeled probe is 
selected to have a first nucleic acid sequence which is 
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complementary to a first portion of a target nucleic acid and 
wherein the nucleic acid sequence of at least one probe 
affixed to the substrate is complementary to a second portion 
of the nucleic acid sequence of the target, the second portion 
being adjacent to the first portion; applying a target nucleic 
acid to the array under suitable conditions for hybridization 
of probe sequences to complementary sequences; introducing a 
labeled probe to the array; hybridizing a probe affixed to the 
substrate to the target nucleic acid; hybridizing the labeled 
probe to the target nucleic acid; affixing the labeled probe to 
an adjacently hybridized probe in the array; and detecting the 
labeled probe affixed to the probe in the array. According to 
preferred methods of the invention the array of probes affixed 
to the substrate comprises a universal set of probes. 
According to other preferred aspects of the invention at least 
two of the probes affixed to the substrate define overlapping 
sequences of the target nucleic acid sequence and more 
preferrably at least two of the labelled probes define 
overlapping sequences of the target nucleic acid sequences. 
Still further, according to another aspect of the invention a 
method is provided for detecting a target nucleic acid of 
known sequence comprising the steps of: contacting a nucleic 
acid sample with a set of immobilized oligonucleotide probes 
attached to a solid substrate under hybridizing conditions 
wherein the immobilized probes are capable of specific 
hybridization with different portions of said target nucleic 
acid sequence; contacting the target nucleic acid with a set 
of labelled oligonucleotide probes in solution under 
hybridizing conditions wherein the labeled probes are capable 
of specific hybridization with different portions of said 
target nucleic acid sequence adjacent to the immobilized 
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probes; covalently joining the immobilized probes to labelled 
probes that are immediately adjacent to the immobilized probe 
on the target sequence (e.g., with ligase) ; removing any 
non-ligated labelled probes; detecting the presence of the 
5 target nucleic acid by detecting the presence of said labelled 
probe attached to the immobilized probes. The invention also 
provides a method of determining expression of a member of a 
set of partially or completely sequenced genes in a cell type, 
a tissue or a tissue mixture comprising the steps of : defining 

10 pairs of fixed and labeled probes specific for the sequenced 
gene; hybridizing unlabeled nucleic acid sample and 
corresponding labeled probes to one or more arrays of fixed 
probes; forming covalent bonds between adjacent hybridized 
labeled and fixed probes; removing unligated probes; and 

15 determining the presence of the sequenced gene by detection of 
labeled probes bound to prespecified locations in the array. 
In a preferred embodiment of this aspect of the invention, the 
target nucleic acid will identif y the pr e sence ,..o, f .. an 
infectious agent T 

20 Further, the present invention provides for an array 

of oligonucleotide probes comprising a nylon membrane; a 
plurality of subarrays of oligonucleotide probes on the nylon 
membrane, the subarrays comprising a plurality of individual 
spots wherein each spot is comprised of a plurality of 

25 oligonucleotide probes of the same sequence; and a plurality 
of hydrophobic barriers located between the subarrays on the 
nylon membrane, whereby the plurality of hyydrophobic barriers 
prevents cross contamination between adjacent subarrays. 

Still further, the present invention provides a 
30 method for sequencing a repetitive sequence, having a first 
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end and a second end, in a target nucleic acid comprising the 
steps of: (a) providing a plurality of spacer oligonucleotides 
of varying lengths wherein the spacer oligonucleotides 
comprise the repetitive sequence; (b) providing a first 
5 oligonucleotide that is known to be adjacent to the first end 
of the repetitive sequence; (c) providing a plurality of 
second oligonucleotides one of which is adjacent to the second 
end of the repetitive sequence, wherein the plurality of 
second oligonucleotides is labeled; (d) hybridizing the first 

10 and the plurality of second oligonucleotides, and one of the 
plurality of spacer oligonucleotides to the target nucleic 
acid ; (e) ligating the hybridized oligonucleotides; (f) 
separating ligated oligonucleotides from unligated 
oligonucleotides; and (g) detecting label in the ligated 

15 oligonucleotides . 

Still further, the present invention provides a 
method for sequencing a branch point sequence, having a first 
end and a second end, in a target nucleic acid comprising the 
steps of: (a) providing a first oligonucleotide that is 

20 complementary to a first portion of the branch point sequence 
wherein the first oligonucleotide extends from the first end 
of the branch point sequence by at least one nucleotide; (b) 
providing a plurality of second oligonucleotides that are 
labeled, and are complementary to a second portion of the 

25 branch point sequence wherein the plurality of second 

oligonucleotides extend from the second end of the branch 
point sequence by at least one nucleotide, and wherein the 
portion of the second oligonucleotides that extend from the 
second end of the branch point " sequence comprise sequences 

30 that are complementary to a plurality of sequences that arise. 
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from the branch point sequence; (c) hybridizing the first 
oligonucleotide, and one of the plurality of second 
oligonucleotides to the target DNA; (d) ligating the 
hybridized oligonucleotides; (e) separating ligated 
oligonucleotides from unligated oligonucleotides; and (f) 
detecting label in the ligated oligonucleotides. 

Still further, the present invention provides a 
method for confirming a sequence by using probes that are 
predicted to be negative for the target nucleic acid. The 
sequence of a target is then confirmed by hybridizing the 
target nucleic acid to the "negative" probes to confirm that 
these probes do not form perfect matches with the target 
nucleic acid. 

Still further, the present invention provides a 
method for analyzing a nucleic acid using oligonucleotide 
probes that are complexed with different labels so that the 
probes may be multiplexed in a hybridization reaction without 
a loss of sequence information (i.e., different probes have 
different labels so that hybridization of the different probes 
to the target can be distinguished) . In a preferred 
embodiment, the labels are radioisotopes, or floursecent 
molecules, or enzymes, or electrophore mass labels. In a more 
preferred embodiment, the differently labeled oligonucleotides 
probes are used in format III SBH, and multiple probes (more 
than two, with one probe being the immobilized probe) are 
ligated together. 

Still further, the present invention provides a 
method for detecting the presence of a target nucleic acid 
having a known sequence when the target is present in very 
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small amounts compared to homologous nucleic acids in a 
sample. In a preferred embodiment, the target nucleic acid is 
an allele present at very low frequency in a sample that has 
nucleic acids from a large number of sources* In an 
5 alternative preferred embodiment, the target nucleic acid has 
a mutated sequence, and is present at very low frequency 
within a sample of nucleic acids. 

Still further, the present invention provides a 
method for confirming the sequence of a target nucleic acid by 

10 using single pass gel sequencing. Primers for single pass gel 
sequencing are derived from the seqeunce obtained by SBH, and 
these primers are used in standard Sanger sequencing reactions 
to provide gel sequence information for the target nucleic 
acid. The sequence obtained by single pass gel sequencing is 

15 then compared to the SBH derived sequence to confirm the 
sequence . 

Still further, the present invention provides a 
method for solving branch points by using single pass gel 
sequencing. Primers for the single pass gel sequencing 

20 reactions are identified from the ends of the Sfs obtained 

after a first round of SBH sequencing, and these primers are 
used in standard Sanger- sequencing reactions to provide gel 
sequencing information through the branch points of the Sfs. 
Sfs are then aligned by comparing the Sanger- sequencing 

25 results through the branch points to the Sfs to identify 
adjoining Sfs. 



Still further, the present invention provides for a 
method of preparing a sample containing target nucleic acids 
by PCR, without purifying the PCR products prior to the SBH 
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reactions. In Format I SBH, crude PCR products are applied to 
a substrate without prior purification, and the substrate may 
be washed prior to introduction of the labeled probes. 

Still further, the present invention provides a 
method and an apparatus for analyzing a target nucleic acid. 
The apparatus comprises two arrays of nucleic acids that are 
mixed together at the desired time. In a preferred 
embodiment, the nucleic acids in one of the arrays are 
labeled. In a more preferred embodiment, a material is 
disposed between the two arrays and this material prevents the 
mixing of nucleic acids in the arrays. When this material is 
removed, or rendered permeable, the nucleic acids in the two 
arrays are mixed together. In an alternative preferred 
embodiment, the nucleic acids in one array are target nucleic 
acids and the nucleic acids in the other are oligonucleotide 
probes. In another preferred embodiment, the nucleic, acids in 
both arrays are oligonucleotide probes. In another preferred 
embodiment, the nucleic acids in one array are oligonucleotide 
probes and target nucleic acids, and nucleic acids in the 
other array are oligonucleotide probes. In another preferred 
embodiment, the nucleic acids in both arrays are 
oligonucleotide probes and target nucleic acids. 

One method of the present invention using the 
apparatus described above comprises the steps of providing an 
array of nucleic acids fixed to a substrate, providing a 
second array of nucleic acids, providing conditions that allow 
the nucleic acids in the second array to come into contact 
with the nucleic acids of the fixed array wherein one of the 
arrays of nucleic acids are target nucleic acids and the other 
array is oligonucleotide probes, and analyzing the 
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hybridization results. In a preferred embodiment, the fixed 
array is target nucleic acid and the second array is lableled 
oligonucleotide probes. In a more preferred embodiment, there 
is a material disposed between the two arrays that prevents 
5 mixing of the nucleic acids until the material is removed or 
rendered permeable to the nucleic acids. 

In a second method of the present invention using 
the apparatus described above comprises the steps of providing 
two arrays of nucleic acid probes, providing conditions that 

10 allow the two arrays of probes, to come into contact with each 
other and a target nucleic acid, ligating together probes that 
are adjacent on the target nucleic acid, and analyzing the 
results. In a preferred embodiment, the probes in one array 
are fixed and the probes in the other array are labeled. In a 

15 more preferred embodiment, there is a material disposed 

between the two arrays that prevents mixing of the probes 
until the material is removed or rendered permeable to the 
probes . 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Format 1 SBH is appropriate for the simultaneous 
analysis of a large set of' samples. Parallel scoring of 
thousands of samples on large arrays may be performed in 
thousands of independent hybridization reactions using small 
25 pieces of membranes. The identification of DNA may involve 

1-20 probes per reaction and the identification of mutations 
may in some cases involve more than 1000 probes specifically 
selected or designed for each sample . For identification of 
the nature of the mutated DNA segments, specific probes may be 
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synthesized or selected for each mutation detected in the 
first round of hybridizations. 

DNA samples may be prepared in small arrays which 
may be separated by appropriate spacers, and which may be 
5 simultaneously tested with probes selected from a set of 
oligonucleotides which may be arrayed in multiwell plates. 
Small arrays may consist of one or more samples. DNA . samples 
in each small array may include mutants or individual samples 
of a sequence. Consecutive small arrays may be organized into 

10 larger arrays. Such larger arrays may include replication of 
the same small array or may include arrays of samples of 
different DNA fragments. A universal set of probes includes 
sufficient probes to analyze a DNA fragment with prespecified 
precision, e.g. with respect to the redundancy of reading each 

15 base pair ("bp"). These sets may include more probes than are 
necessary for one specific fragment, but may include fewer 
probes than are necessary for testing thousands of DNA samples 
of different sequence. 

DNA or allele identification and a diagnostic 
sequencing process may include the steps of: 

1) Selection of a subset of probes from a dedicated, 
representative or universal set to be hybridized with 
each of a plurality of small arrays; 

2) Adding a first probe to each subarray on each of the 
arrays to be analyzed in parallel; 

3) Performing hybridization and scoring of the 
hybridization results; 

4) Stripping off previously used probes ; 

5) Repeating hybridizaton, scoring and stripping steps 
for the remaining probes which are to be scored; 



20 



25 



30 
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5) 



Processing the obtained results to obtain a final 



analysis or to determine additional probes to be 



hybridized; 



6) 



Performing additional hybridizations for certain 



5 



subarrays ; and 



7) Processing complete sets of data and obtaining a 
final analysis. 

This approach provides fast identification and 
sequencing of a small number of nucleic acid samples of one 

10 tv P e (e.g. DNA, RNA) , and also provides parallel analysis of 
many sample types in the form of subarrays by using a 
presynthesized set of probes of manageable size. Two 
approaches have been combined to produce an efficient and 
versatile process for the determination of DNA identity, for 

15 DNA diagnostics, and for identification of mutations. 



set of shorter probes may be used in place of a longer unique 
probe. In this approach, although there may be more probes to 
be scored, a universal set of probes may be synthesized to 
20 cover any type of sequence. For example, a full set of 6-mers 
includes only 4,096 probes, and a complete set of 7-mers 
includes only 16,384 probes. 



with two levels of hybridization. One level is hybridization 
25 of a sufficient set of probes that cover every base at least 
once. For this purpose, a specific set of probes may be 
synthesized for a standard sample. The results of 
hybridization with such a set of probes reveal whether and 
where mutations (differences) occur in non-standard samples. 



For the identification of known sequences, a small 



Full sequencing of a DNA fragment may be performed 



WO 98/31836 PCTYUS98/00704 

13 

Further, this set of probes may include "negative" probes to 
confirm the hybridization results of the "positive" probes. 
To determine the identity of the changes, additional specific 
probes may be hybridized to the sample . This additional set 
5 of probes will have both "positive" (the mutant sequence) and 
"negative" probes, and the sequence changes will be identified 
by the positive probes and confirmed by the negative probes. 

In another embodiment, all probes from a universal 
set may be scored. A universal set of probes allows scoring 
of a relatively small number of probes per sample in a two 
step process without an undesirable expenditure of time. The 
hybridization process may involve successive probings, in a 
first step of computing an optimal subset of probes to be 
hybridized first and, then, on the basis of the obtained 
results, a second step of determining additional probes to be 
scored from among those in a universal set. Both sets of 
probes have "negative" probes that confirm the positive probes 
in the set. Further, the sequence that is obtained may then 
be confirmed in a spearate step by hybridizing the sample with 
a set of "negative" probes identified from the SBH results. 

In SBH sequence assembly, K -1 oligonucleotides 
which occur repeatedly in analyzed DNA fragments due to chance 
or biological reasons may be subject to special consideration. 
If there is no additional information, relatively small 
25 fragments of DNA may be fully assembled in as much as every 
base pair is read several times* 

In the assembly of relatively longer fragments, 
ambiguities may arise due to the repeated occurrence in a set 
of positively-scored probes of a K-l sequence (i.e., a 



15 
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sequence shorter than the length of the probe) . This problem 
does not exist if mutated or similar sequences have to be 
determined (i.e., the K-l sequence is not identically 
repeated) . Knowledge of one sequence may be used as a template 
to correctly assemble a sequence known to be similar (e.g. by 
its presence in a database) by arraying the positive probes 
for the unknown sequence to display the best fit on the 
template . 

The use of an array of sample avoids consecutive 
scoring of many oligonucleotides on a single sample or on a 
small set of samples. This approach allows the scoring of more 
probes in parallel by manipulation of only one physical 
object. Subarrays of DNA samples 1000 bp in length may be 
sequenced in a relatively short period of time. If the samples 
are spotted at 50 subarrays in an array and the array is 
reprobed 10 times, 50 0 probes may be scored. In screening for 
the occurrence of a mutation, enough probes may be used to 
cover each base three times. If a mutation is present, several 
covering probes will be affected. The use of information about 
the identity of negative probes may map the mutation with a 
two base precision. To solve a single base mutation mapped in 
this way, an additional 15 probes may be employed. These 
probes cover any base combination for two questionable 
positions (assuming that deletions and insertions are not 
involved) . These probes may be scored in one cycle on 5 0 
subarrays which contain a given sample. In the implementation 
of a multiple label color scheme (i.e., multiplexing), two to 
six probes, each having a different label such as a different 
fluorescent dye, may be used as a pool, thereby reducing the 
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number of hybridization cycles and shortening the sequencing 
process . 

In more complicated cases, there may be two close 
mutations or insertions. They may be handled with more probes. 
For example, a three base insertion may be solved with 64 
probes . The most complicated cases may be approached by 
several steps of hybridization, and the selecting of a new set 
of probes on the basis of results of previous hybridizations. 

If subarrays to be analyzed include tens or hundreds 
of samples of one type, then several of them may be found to 
contain one or more changes (mutations, insertions, or 
deletions) . For each segment where mutation occurs, a specific 
set of probes may be scored. The total number of probes to be 
scored for a type of sample may be several hundreds . The 
scoring of replica arrays in parallel facilitates scoring of 
hundreds of probes in a relatively small number of cycles. In 
addition, compatible probes may be pooled. Positive 
hybridizations may be assigned to the probes selected to check 
particular DNA segments because these segments usually differ 
in 75% of their constituent bases. 

By using a larger set of longer probes, longer 
targets may be analyzed. These targets may represent pools of 
fragments such as pools of exon clones. 

A specific hybridization scoring method may be 
employed to define the presence of mutants in a genomic 
segment to be sequenced from a diploid chromosomal set. Two 
variations are where: i) the sequence from one chromosome 
represents a known allele and the sequence from the other 
represents a new mutant; or, ii) both chromosomes contain new, 
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but different mutants. In both cases, the scanning step 
designed to map changes gives a maximal signal difference of 
two -fold at the mutant position. Further, the method can be 
used to identify which alleles of a gene are carried by an 
5 individual and whether theindividual is homozygous or 
heterozygous for that gene . 

Scoring two -fold signal differences required in the 
first case may be achieved efficiently by comparing 
corresponding signals with homozygous and heterozygous 

10 controls. This approach allows determination of a relative 
reduction in the hybridization signal for each particular 
probe in a given sample. This is significant because 
hybridization efficiency may vary more than two- fold for a 
particular probe hybridized with different nucleic acid 

15 fragments having the same full match target. In addition, 
different mutant sites may affect more than one probe 
depending upon the number of oligonucleotide probes. Decrease 
of the signal for two to four consecutive probes produces a 
more significant indication of a mutant site. Results may be 

20 checked by testing with small sets of selected probes among 
which one or few probes selected to give a full match signal 
which is on average eight -fold stronger than the signals 
coming from mismatch- containing duplexes. 

Partitioned membranes allow a very flexible 
25 organization of experiments to accommodate relatively larger 

numbers of samples representing a given sequence type, or many 
different types of samples represented with relatively small 
numbers of samples. A range of 4-256 samples can be handled 
with particular efficiency. Subarrays within this range of 
30 , numbers of dots may be designed to match the configuration and 
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size of standard multiwell plates used for storing and 
labelling oligonucleotides. The size of the subarrays may be 
adjusted for different number of samples, or a few standard 
subarray sizes may be used. If all samples of a type do not 
5 fit in one subarray, additional subarrays or membranes may be 
used and processed with the same probes. In addition, by 
adjusting the number of replicas for each subarray, the time 
for completion of identification or sequencing process may be 
varied . 

10 As used herein, "intermediate fragment" means an 

oligonucleotide between 5 and 1000 bases in length, and 
preferably between 10 and 40 bp in length. 

In Format 3, a first set of oligonucleotide probes 
of known sequence is immobilized on a solid support under 

15 conditions which permit them to hybridize with nucleic acids 
having respectively complementary sequences. A labeled, 
second set of oligonucleotide probes is provided in solution. 
Both within the sets and between the sets the probes may be 
of the same length or of different lengths. A nucleic acid to 

20 be sequenced or intermediate fragments thereof may be applied 
to the first set of probes in double -stranded form 
(especially where a recA protein is present to permit 
hybridization under non- denaturing conditions) , or in 
single -stranded form and under conditions which permit hybrids 

25 of different degrees of complementarity (for example, under 

conditions which allow discriminatation between full match and 
one base pair mismatch hybrids) . The nucleic acid to be 
sequenced or intermediate fragments thereof may be applied to 
the first set of probes before, after or simultaneously with 

30 the second set of probes. Probes that bind to adjacent sites 
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on the atrget are bound together (e.g., by stacking 
interactions or by a ligase or other means of causing 
chemical bond formation between the adjacent probes) . After 
permitting adjacent probes to be bound, fragments and probes 
5 which are not immobilized to the surface by chemical bonding 
to a member of the first set of probe are washed away, for 
example, using a high temperature (up to 100 degrees C) wash 
solution which melts hybrids. The bound probes from the 
second set may then be detected using means appropriate to the 
10 label employed (which may, for example, be chemiluminescent , 
fluorescent, radioactive, enzymatic, densitometric , or 
electrophore mass labels) . 

Herein, nucleotide bases "match" or are 
"complementary" if they form a stable duplex by hydrogen 

15 bonding under specified conditions. For example, under 

conditions commonly employed in hybridization assays, adenine 
("A") matches thymine ( M T ,! ), but not guanine ("G") or cytosine 
("C") . Similarly, G matches C, but not A or T\ Other bases 
which will hydrogen bond in less specific fashion, such as 

20 inosine or the Universal Base ("M" base, Nichols et al 1994), 
or other modified bases, such as methylated bases, for 
example, are complementary to those bases for which they form 
a stable duplex under specified conditions. A probe is said 
to be "Perfectly complementary" or is said to be a "p e-r-F^nt- 

25 match" if , 5ach b ase in the prob e forms a duplex by hydrogen 
bonding to a base in the nucleic an I d to be sequenced 
according to the Watson and Crick b a se paring rule.cs fj. fi . r 
absent any surroun ding frequence e ffects, the duplex -Formed has 
the maximal binding energy for a p arti cular probe). 

30 "Perfectly complementary" and "perfect match" a re a i ar > 
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to encompass probes which have analogs or modified 
n ucleotides. A "perfect match" for an analog or modified 
nucleotide is judged according to a "perfect match rule" 
selected for that analog or modified nucleotide (e.g., the 
5 binding pair that has maximal binding energy for a particular 

analog or modified nucle otide) . Bach base in a probe that 

does not form a binding pair according to the "rules" is said 
to be a "mismatch" under the specified hybridization 
conditions . 

10 A list of probes may be assembled wherein each probe 

is a perfect match to the nucleic acid to be sequenced. The 
probes on this list may then be analyzed to order them in 
maximal overlap fashion. Such ordering may be accomplished by 
comparing a first probe to each of the other probes on the 

15 list to determine which probe has a 3 ' end which has the 

longest sequence of bases identical to the sequence of bases 
at the 5' end of a second probe. The first and second probes 
may then be overlapped; and the process may be repeated by 
comparing the 5' end of the second probe to the 3' end of all 

20 of the remaining probes and by comparing the 3 1 end of the 

first probe with the 5' end of all of the remaining probes. 
The process may be continued until there are no probes on the 
list which have not been overlapped with other probes. 
Alternatively, more than one probe may be selected from the 

25 list of positive probes, and more than one set of overlapped 

probes ("sequence nucleus") may be generated in parallel. The 
list of probes for either such process of sequence assembly 
may be the list of all probes which are perfectly 
complementary to the nucleic acid to be sequenced or may be 

30 any subset thereof . 
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The 5' and 3' ends of the probes may be overlapped 
to generate longer stretches of sequence. This process of 
assembling probes continues until an ambiguity arises because 
of a branch point (a probe is repeated in the fragment) , 
5 repetitive sequences longer than the probes, or an uncloned 

segment. The streches of sequence between any two ambiguities 
are referred to as fragment os a subclone sequence (Sfs) . 
Where ambiguities arise in sequence assembly due to the 
availability of alternative proper overlaps with probes, 
10 hybridization with longer probes spanning the site of overlap 
alternatives, competitive hybridization, ligation of 
alternative end to end pairs of probes spanning the site of 
ambiguity or single pass gel analysis (to provide an 
unambiguous ordering of Sfs) may be used. 

15 By employing the above procedures, one may obtain 

any desired level of sequence, from a pattern of 
hybridization (which may be correlated with the identity of a 
nucleic acid sample to serve as a signature for identifying 
the nucleic acid sample) to overlapping or non- over lapping 

20 probes up through assembled Sfs and on to complete sequence 

for an intermediate fragment or an entire source DNA molecule 
(e.g. a chromosome) . 

Sequencing may generally comprise the following 

steps : 

25 (a) contacting an array of immobilized 

oligonucleotide probes with a nucleic acid fragment under 
conditions effective to allow the fragment to form a primary 
complex with an immobilized probe having a complementary 
sequence ; 
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(b) contacting this primary complex with a set of 
labeled oligonucleotide probes in solution under conditions 
effective to allow the primary complex to hybridize to the 
labeled probe, thereby forming secondary complexes wherein the 

5 fragment is hybridized with both an immobilized probe and a 
labeled probe; 

(c) removing from a secondary complex any labeled 
probe that has not hybridized adjacent to an immobilized 
probe ; 

10 (d) detecting the presence of adjacent labeled and 

unlabeled probes by detecting the presence of the label; and 

(e) determining a nucleotide sequence of trie 
fragment by connecting the known sequence of the immobilized 
and labeled probes. 

15 Hybridization and washing conditions may be selected 

to detect substantially perfect match hybrids (such as those 
wherein the fragment and probe hybridize at six out of seven 
positions) / may be selected to allow differentiation of 
perfect matches and one base pair mismatches, or may be 

20 selected to permit detection only of perfect match hybrids. 

Suitable hybridization conditions may be routinely 
determined by optimization procedures or pilot studies* Such 
procedures and studies are routinely conducted by those 
skilled in the art to establish protocols for use in a 
25 laboratory. See e.g., Ausubel et al . , Current Protocols in 
Molecular Biology, Vol, 1-2, John Wiley & Sons (1989) ; 
Sambrook et al . , Molecular Cloning A Laboratory Manual, 2nd 
Ed., Vols. 1-3, Gold Springs Harbor Press (1989); and Maniatis 
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et al., Molecular Cloning: A Laboratory Manual , Cold Spring 
Harbor Laboratory Cold Spring Harbor, New York (1982) , all of 
which are incorporated by reference herein. For example, 
conditions such as temperature, concentration of components, 
5 hybridization and washing times, buffer components, and their 
pH and ionic strength may be varied. 

In embodiments wherein the labeled and immobilized 
probes are not physically or chemically linked, detection may 
rely solely on washing steps of controlled stringency. Under 
10 such conditions, adjacent probes have increased binding 

affinity because of stacking interactions between the adjacent 
probes. Conditions may be varied to optimize the process as 
described above . 

In embodiments wherein the immobilized and labeled 
probes are ligated, ligation may be implemented by a chemical 
ligating agent {e.g. water-soluble carbodiimide or cyanogen 
bromide) , or a ligase enzyme, such as the commercially 
available T 4 DNA ligase may be employed. The washing 
conditions may be selected to distinguish between adjacent 
versus nonadjacent labeled and immobilized probes exploiting 
the difference in stability for adjacent probes versus 
nonadjacent probes. 

Oligonucleotide probes may be labeled with 
fluorescent dyes, chemiluminescent systems, radioactive labels 
25 (e.g., 35 S, 3 H, 32 P or 33 P) or with isotopes detectable by mass 
spectrometry . 

Where a nucleic acid molecule of unknown sequence is 
longer than about 4 5 or 50 bp, the molecule may be fragmented 
and the sequences of the fragments determined. Fragmentation 
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may be accomplished by restriction enzyme digestion, shearing 
or NaOH. Fragments may be separated by size (e.g. by gel 
electrophoresis) to obtain a preferred fragment length of 
about ten to forty bps . 

5 Oligonucleotides may be immobilized, by a number of 

methods known to those skilled in the art, such as 
laser-activated photodeprotection attachment through a 
phosphate group using reagents such as a nucleoside 
phosphoramidite or a nucleoside hydrogen phosphorate. Glass, 

10 nylon, silicon and fluorocarbon supports may be used. 

Oligonucleotides may be organized into arrays, and 
these arrays may include all or a subset of all probes of a 
given length, or sets of probes of selected lengths. 
Hydrophobic partitions may be used to separate probes or 

15 subarrays of probes. Arrays may be designed for various 

applications (e.g. mapping, partial sequencing, sequencing of 
targeted regions for diagnostic purposes, mRNA sequencing and 
large scale sequencing) . A specific chip may be designed to be 
dedicated to a particular application by selecting a 

20 combination and arrangement of probes on a substrate* 

For example, 1024 immobilized probe arrays of all 
oligonucleotide probes 5 bases in length (each array 
containing 10 24 distinct probes) may be constructed. The 
probes in this example are 5-mers in an informational sense 
25 (they may actually be longer probes) ♦ A second set of 1024 5- 
mer probes may be labeled, and one of each labeled probe may 
be applied to an array of immobilized probes along with a 
fragment to be sequenced. In this example, 1024 arrays would 
be combined in a large superarray, or "superchip . " In those 
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instances where an immobilized probe and one of the labeled 
probes hybridize end -to-end along a nucleic acid fragment., 
the two probes are joined, for example by ligation, and, after 
removing unbound label, 10-mers complementary to the sample 
5 fragment are detected by the correlation of the presence of a 
label at a point in an array having an immobilized probe .of 
known sequence to which was applied a labeled probe of known 
sequence. The sequence of the sample fragment is simply the 
sequence of the immobilized probe continued in the sequence 
10 of the labeled probe. In this way, all one million possible 
10-mers may be tested by a combinatorial process which employs 
only 5-mers and which thus involves one thousandth of the 
amount of effort for oligonucleotide synthesis. 

A nucleic acid sample to be sequenced may be 
fragmented or otherwise treated (for example, by the use of 
recA) to avoid hindrance. to hybridization from secondary 
structure in the sample. The sample may be fragmented by, for 
example, digestion with a restriction enzyme such as Cvi ji , 
physical shearing (e.g. by ultrasound ), or by NaOH treatment. 
The resulting fragments may be separated by gel 
electrophoresis and fragments of an appropriate length, such 
as between about 10 bp and about 4 0 bp, may be extracted from 
the gel . 

A reusable Format 3 SBH array may be produced by 
introducing a cleavable bond between the fixed and labeled 
probes and then cleaving this bond after a round of Format 3 
analyziz is finished. The labeled probes may be 
ribonucleotides or a ribonucleotide may be used as the joining 
base in the labeled probe so that this probe may subsequently 
be removed, e.g., by RNAse'or uracil-DNA glycosylate 
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treatment, or NaOH treatment. In addition, bonds produced by 
chemical ligation may be selectively cleaved. 

Other variations include the use of modified 
oligonucleotides to increase specificity or efficiency, 

5 cycling hybridizations to increase the hybridization signal, 
for example by performing a hybridization cycle under 
conditions (e.g. temperature) optimally selected for a first 
set of labeled probes followed by hybridization under 
conditions optimally selected for a second set of labeled 

10 probes. Shifts in reading frame may be determined by using 

mixtures (preferably mixtures of eguimolar amounts) of probes 
ending in each of the four nucleotide bases A, T, C and G. 

Branch points produce ambiguities as to the ordered 
sequence of a fragment . Although the sequence information is 

15 determined by SBH, either: (i) long read length, single-pass 
gel sequencing at a fraction of the cost of complete gel 
sequencing; or (ii) comparison to related sequences, may be 
used to order hybridization data where such ambiguities 
{ "branch points' 7 ) occur. Primers for single pass gel 

20 sequencing throug h the branch poin ts ar e identified from .. the . 
SRH sequenc e information or. from known vector seguencep, e,g,, 
the flanking geguen c eK to the vector insert site, and Standard 
flan gfir'-sfiquenflina re actions are performed on the sample 
nucleic acid. The s equence obtained from this single pass gel 

25 sequencing is com pared to the Sfs that read into and QUt of 
the branch p oints to identify the order of the SfS, 
Alternatively, the S fs may be ordered by comparing the 
sequence of the Sffi to related se quences and ordering the Sfs 
to produce a sequenc e that is closest to the related sequence,,,,, 
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In addition, the number of tandem repetitive nucleic 
acid segments in a target fragment may be determined by 
single-pass gel sequencing. As tandem repeats occur rarely in 
protein-encoding portions of a gene, the gel -sequencing step 
5 will be performed only when one of these noncoding regions is 
identified as being of particular interest (e.g., if it is an 
important regulatory region) . 

Obtaining information about the degree of 
hybridization exhibited for a set of only about 200 

10 oligonucleotides probes (about 5% of the effort required for 
complete sequencing) defines a unique signature of each gene 
and may be used for sorting the cDNAs from a library to 
determine if the library contains multiple copies of the same 
gene. By such signatures, identical, similar and different 

15 cDNAs can be distinguished and inventoried. 

Nucleic acids and methods for isolating, cloning and 
sequencing nucleic acids are well known to those of skill in 
the art. See e.g., Ausubel et al., Current Protocols in 
Molecular Biology, Vol. 1-2, John Wiley & Sons (1989); and 
20 Sambrook et al., Molecular Cloning A Laboratory Manual, 2nd 
Ed., Vols. 1-3, Cold Springs Harbor Press (1989), both of 
which are incorporated by reference herein. 

SBH is a well developed technology that may be 
practiced by a number of methods known to those skilled in the 
25 art. Specifically, techniques related to sequencing by 

hybridization of the following documents is incorporated by 
reference herein: Drmanac et al . , U.S. Patent No. 5,202,231 
(hereby incorporated by reference herein) - Issued April 13, 
1993; Drmanac et al . , Genomics, 4, 114-128 (1989); Drmanac et 
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al . , Proceedings of the First Int ' 1 . Conf. Electrophoresis 
Super computing Human Genome Cantor et al . eels, World 
Scientific Pub. Co., Singapore, 47-59 (1991); Drmanac et al . , 
Science, 260, 1649-1652 (1993); Lehrach et al . , Genome 
5 Analysis; Genetic and Physical Mapping, 1, 39-81 (1990) , Cold 
Spring Harbor Laboratory Press; Drmanac et al . , Nucl . Acids 
Res., 4691 (1986); Stevanovic et al . , Gene, 79, 139 (1989); 
Panusku et al . , Mol . Biol. Evol . , 1, 607 (1990); Nizetic et 
al . , Nucl. Acids Res . t 19, 182 (1991); Drmanac et al . , J. 
10 Biomol. Struct. Dyn. , 5, 1085 (1991); Hoheisel et al . , Mol. 

Gen., 4, 125-132 (1991); Strezoska et al . , Proc . Nat ' 1 . Acad. 
Sci. (USA), 88, 10089 (1991); Drmanac et al . , Nucl. Acids 
Res., 19, 5839 (1991); and Drmanac et al . , Jnt. J. Genome 
Res. , 1, 59-79 (1992) . 

15 

The present invention is illustrated in the 
following examples. Upon consideration of the present 
disclosure, one of skill in the art will appreciate that many 
other embodiments and variations may be made in the scope of 
20 the present invention. Accordingly, it is intended that the 

broader aspects of the present invention not be limited to the 
disclosure of the following examples. 
EX AM PLE I 

Preparation of Sets of Probes 

25 Two types of universal sets of probes may be 

prepared. The first is a complete set (or at least a 
noncomplementary subset) of relatively short probes, for 
example all 4096 (or about 2000 non- complementary) 6-mers, or 
all 16,384 (or about 8,000 non-complementary) 7-mers. Full 
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noncomplementary subsets of 8-mers and longer probes are less 
convenient inasmuch as they include 32,000 or more probes, 

A second type of probe set is selected as a small 
subset of probes still sufficient for reading every bp in any 
5 sequence with at least with one probe. For example, 12 of 16 
dimers are sufficient. A small subset for 7-mers, 8-mer and 
9-mers for sequencing double stranded DNA may be about 3000, 
10,0 00 and 3 0,00 0 probes, respectively. 

Sets of probes may also be selected to identify a 
10 target n ucleic acid of known sequence, and/or to identify 
al ifilea or mutants of a target nucleic acid with a known 
sequence. Such a set of probes contains sufficient probes so 
that every nucleotide position of the target nucleic acid is 

read at least once, Al l eles .,. Qr m utant s are identified by the 

15 loss of binding of one of the "positive" probes. The specific 

sequence of these alleles or mutants is then determined bv 
interrogating the target nucleic acid wit h sets of probes that 
contain every possible nucleotide chang e and combination of 
Changes at these* probe positions t 

Probes may be prepared using standard chemistry with 
one to three non-specified (mixed A,T,C and G) or universal 
(e.g. M base or inosine) bases at the ends. If radiolabelling 
is used, probes may have an OH group at the 5 » end for 
kinasing by radiolabelled phosphorous groups. Alternatively, 
probes labelled with any compatible system, such as 
fluorescent dyes, may be employed. Other types of probes, such 
as PNA (Protein Nucleic Acids) or probes containing modified 
bases which change duplex stability also may be used. 
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Probes may be stored in bar-coded multiwell plates. 
For small numbers of probes , 96 -well plates may be used; for 
10,000 or more probes, storage in 384- or 864-well plates is 
preferred. Stacks of 5 to 5 0 plates are enough to store all 
5 probes. Approximately 5 pg of a probe may be sufficient for 
hybridization with one DNA sample. Thus, from a small 
synthesis of about 50 mg per probe, ten million samples may be 
analyzed. If each probe is used for every third sample, and if 
each sample is 1000 bp in length, then over 30 billion bases 
10 (10 human genomes) may be sequenced by a set of 5,000 probes. 

Probes Having Modified Oliaomicleotideg 

Modified oligonucleotides may be introduced into 
hybridization probes and used under appropriate conditions 

15 therefor. For example, pyrimidines with a halogen at the 
exposition may be used to improve duplex stability by 
influencing base stacking- 2 , 6 -diaminopurine may be used to 
provide a third hydrogen bond in base pairing with thymine, 
thereby thermally stabilizing DNA-duplexes . Using 

20 2 , 5 -diaminopurine may increase duplex stability to allow more 
stringent conditions for annealing, thereby improving the 
specificity of duplex formation, suppressing background 
problems and permitting the use of shorter oligomers. 

The synthesis of the triphosphate versions of these 
25 modified nucleotides is disclosed by Hoheisel & Lehrach 
(1990) . 

One may also use the non-discriminatory base 
analogue, or universal base, as designed by Nichols et al ♦ 
(19 94) ♦ This new analogue, l- (2 -deoxy- 
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-D-ribf uranosyl) ~3 -nitropyrrole (designated M) , was generated 
for use in oligonucleotide probes and primers for solving the 
design problems that arise as a result of the degeneracy of 
the genetic code, or when only fragmentary peptide sequence 
5 data are available. This analogue maximizes stacking while 
minimizing hydrogen -bonding interactions without sterically 
disrupting a DNA duplex. 

The M nucleoside analogue was designed to maximize 
stacking interactions using aprotic polar substituents linked 
to heteroaromatic rings, enhancing intra- and inter-strand 
stacking interactions to lessen the role of hydrogen bonding 
in base-pairing specificity. Nichols et al . (1994) favored 
3 -nitropyrrole 2 -deoxyribonucleoside because of its 
structural and electronic resemblance to p-nitroaniline , whose 
derivatives are among the smallest known intercaltors of 
double -stranded DNA. 

The dimethoxytrityl-protected phosphoramidite of 
nucleoside M is also available for incorporation into 
nucleotides used as primers for sequencing and polymerase 
20 chain reaction (PCR) . Nichols et al . (1994) showed that a 

substantial number of nucleotides can be replaced by M without 
loss of primer specificity. 

A unique property of M is its ability to replace 
long strings of contiguous nucleosides and still yield 
25 functional sequencing primers. Sequences with three, six and 
nine M substitutions have all been reported to give readable 
sequencing ladders, and PCR with three different M-containing 
primers all resulted in amplification of the correct product 
(Nichols et al . , 1994). 
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The ability of 3 -nitropyrrole-containing 
oligonucleotides to function as primers strongly suggests that 
a duplex structure must form with complementary strands. 
Optical thermal profiles obtained for the oligonucleotide 
pairs d(5 -C 2 -T 5 XT- 5 G 2 -3 ) and d(5 -C 2 A 5 YA 5 G2 -3 ) (where X and Y 
can be A, C, G, T or M) were reported to fit the normal 
sigmoidal pattern observed for the DNA double-to single strand 
transition. The Tm values of the oligonucleotides containing 
X M base pairs (where X was A, C, G or T, and Y was M) were 
reported to all fall within a 3°C range (Nichols et al . , 1994) . 

Selection and Labeling of Probes 

When an array of subarrays is produced, the sets of 
probes to be hybridized in each of the hybridization cycles on 
each of the subarrays is defined. For example, a set of 3 84 
probes may be selected from the universal set, and 96 probings 
may be performed in each of 4 cycles. Probes selected to be 
hybridized in one cycle preferably have similar G+C contents. 

Selected probes for each cycle are transferred to a 
96 -well plate and then are labelled by kinasing or by other 
labelling procedures if they are not labelled (e.g. with 
stable fluorescent dyes) before they are stored. 

On the basis of the first round of hybridizations, a 
new set of probes may be defined for each of the subarrays for 
additional cycles. Some of the arrays may not be used in some 
of the cycles. For example, if only 8 of 64 patient samples 
exhibit a mutation and 8 probes are scored first for each 
mutation, then all 64 probes may be scored in one cycle and 3 2 
subarrays are not used. These unused subarrays may then be 
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treated with hybridization buffer to prevent drying of the 
filters . 

Probes may be retrieved from the storing plates by 
any convenient approach, such as a single channel pipetting 
5 device, or a robotic station, such as a Beckman Biomek 1000 
(Beckman Instruments, Fullerton, California) or a Mega Two 
robot (Megamation, Lawrenceville, New Jersey) . A robotic 
station may be integrated with data analysis programs and 
probe managing programs. Outputs of these programs may be 
10 inputs for one or more robotic stations. 

Probes may be retrieved one by one and added to 
subarrays covered by hybridization buffer. It is preferred 
that retrieved probes be placed in a new plate and labelled or 
mixed with hybridization buffer. The preferred method of 

15 retrieval is by accessing stored plates one by one and 

pipetting (or transferring by metal pins) a sufficient amount 
of each selected probe from each plate to specific wells in an 
intermediary plate. An array of individually addressable 
pipettes or pins may be used to speed up the retrieval 

20 process. 
ramPM ± 

Preparation of frabo led l Probes 

The oligonucleotide probes may be prepared by 
automated synthesis, which is routine to those of skill in the 
25 art, for example, using and Applied Biosystems system. 
Alternatively, probes may be prepared using Genosys 
Biotechnologies Inc. Methods using stacks of porous Teflon 
wafers . 
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Oligonucleotide probes may be labeled with, for 

3 c 32 33 33 

example, radioactive labels ( S, P, P, and preferably, P) 
for arrays with 100-200 um or 100-400 urn spots; 
non- radioactive isotopes (Jacobsen et al . , 1990); or 
5 fluorophores (Brumbaugh et a.1 . , 1988). All such labeling 

methods are routine in the art, as exemplified by the relevant 
sections in Sambrook et al . (1989) and by further references 
such as Schubert et al . (1990), Murakami et al . (1991) and 
Gate et al . (1991), all articles being specifically 
10 incorporated herein by reference. 

In regard to radiolabelling, the common methods are 
end-labeling using T4 polynucleotide kinase or high specific 
activity labeling using Klenow or even T7 polymerase. These 
are described as follows . 

15 Synthetic oligonucleotides are synthesized without a 

phosphate group at their 5 termini and are therefore easily 
labeled by transfer of the - 32 P or - 33 P from [ - 32 P]ATP or [ 
- 33 P]ATP using the enzyme bacteriophage T4 polynucleotide 
kinase. If the reaction is carried out efficiently, the 

20 specificity activity of such probes can be as high as the 

specific activity of the [ - 32 P]ATP or [ - 33 P]ATP itself. The 
reaction described below is designed to label 10 pmoles of an 
oligonucleotide to high specific activity. Labeling of 
different amounts of oligonucleotide can easily be achieved by 

25 increasing or decreasing the size of the reaction, keeping the 
concentrations of all components constant. 



A reaction mixture would be created using 1.0 ul of 
oligonucleotide (10 pmoles /ul ) ; 2 . 0 ul of 10 x bacteriophage 
T4 polynucleotide kinase buffer; 5.0 ul of [ - 32 P] ATP or [ 
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- 33 P]ATP (sp. Act. 5000 Ci/mmole; 10 mCi/ml in aqueous 
solution) (10 pmoles) ; and 11.4 ul of water* Eight (8) units 
{~1 ul) of bacteriophage T4 polynucleotide kinase is added to 
the reaction mixture, and incubated for 4 5 minutes at 3 7°C. 
5 The reaction is heated for 10 minutes at 68°C to inactivate the 
bacteriophage T4 polynucleotide kinase. 

The efficiency of transfer of 32 P or 33 P to the 
oligonucleotide and its specific activity is then determined. 
If the specific activity of the probe is acceptable, it is 
10 purified. If the specific activity is too low, an additional 
8 units of enzyme is added and incubated for a further 3 0 
minutes at 3 7°C before heating the reaction for 10 minutes at 
6Q°C to inactivate the enzyme. 

Purification of radiolabeled oligonucleotides can be 
15 achieved by, e.g., precipitation with ethanol; precipitation 
with cetylpyridinium bromide; by chromatography through 
bio-gel P-60; or by chromatography on a Sep-Pak C 18 column, or 
by polyacryiamide gel electrophoresis. 

Probes of higher specific activities can be obtained 
20 using the Klenow fragment of E. coli. DNA polymerase I to 
synthesize a strand of DNA complementary to the synthetic 
oligonucleotide. A short primer is hybridized to an 
oligonucleotide template whose sequence is the complement of 
the desired radiolabeled probe. The primer is then extended 
25 using the Klenow fragment of E . coli DNA polymerase I to 
incorporate [ - 32 P] dNTPs or [ - 33 p] dNTPs in a 
template-directed manner. After the reaction, the template 
and product are separated by denaturation followed by 
electrophoresis through a polyacryiamide gel under denaturing 
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conditions. With this method, it is possible to generate 
oligonucleotide probes that contain several radioactive atoms 
per molecule of oligonucleotide. 

To use this method, one would mix in a microfuge 
tube the calculated amounts of [a-32P]dNTPs or [a-33P]dNTPs 
necessary to achieve the desired specific activity and 
sufficient to allow complete synthesis of all template 
strands. Then add to the tube the appropriate amounts of 
primer and template DNAs , with the primer being in three- to 
tenfold molar excess over the template. 

0.1 volume of 10 x Klenow buffer would then be added 
and mixed well. 2-4 units of the Klenow fragment of E.aoll 
DNA polymerase I would then be added per 5 ul of reaction 
volume, mixed and incubated for 2-3 hours at 4°C. If desired, 
the process of the reaction may be monitored by removing small 
(0,1 ul) aliquots and measuring the proportion of 
radioactivity that has become precipitable with 10% 
trichloroacetic acid (TCA) . 

The reaction would be diluted with an equal volume 
of gel-loading buffer, heated to 80oC for 3 minutes, and then 
the entire sample loaded on a denaturing polyacrylamide gel. 
Following electrophoresis, the gel is autoradiographed, 
allowing the probe to be localized and removed from the gel. 
Various methods for fluorescent probe labeling are also 
available, e.g., Brumbaugh et al . (1988) describe the 
synthesis of f luorescently labeled primers. A deoxyuridine 
analog with a primary amine "linker arm" of 12 atoms attached 
at C-5 is synthesized. Synthesis of the analog consists of 
derivatizing 2 -deoxyuridine through organometallic 
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intermediates to give 5 (methyl propenoyl) -2 -deoxyuridine . 
Reaction with dimethoxytrityl- chloride' produces the 
corresponding 5 -dimethoxytrityl adduct. The methyl ester is 
hydrolyzed, activated, and reacted with an appropriately 
5 monoacylated alkyl diamine. After purification, the resultant 
linker arm nucleosides are converted to nucleoside analogs 
suitable for chemical oligonucleotide synthesis. 

Oligonucleotides would then be made that include one 
or two linker arm bases by using modified phosphoridite 

10 chemistry. To a solution of 50 nmol of the linker arm 

oligonucleotide in 25 ul of 500 mM sodium biocarbonate (pH 
9.4) is added 20 ul of 300 mM FITC in dimethyl sulfoxide. The 
mixture is agitated at room temperature for 6 hrs . The 
oligonucleotide is separated from free FITC by elution form a 

15 1 x 30 cm Sephadex G-25 column with 20 mM ammonium acetate (pH 
6) , combining fractions in the first UV-absorbing peak. 

In general, fluorescent labeling of an 
oligonucleotide at its 5 1 -end initially involved two steps . 
First, a iV-protected aminoalkyl phosphoramidite derivative is 

20 added to the 5 ' -end of an oligonucleotide during automated 
nucleic acid synthesis. After removal of all protecting 
groups, the NHS ester of an appropriate fluorescent dye is 
coupled to the 5' -amino group overnight followed by 
purification of the labeled oligonucleotide from the excess of 

25 dye using reverse phase HPLC or PAGE. 



Schubert et al . (1990) described the synthesis of a 
phosphoramidite that enables oligonucleotides labeled with 
fluorescein to be produced during automated DNA synthesis. 
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Murakami et a.1. also described the preparation of 
flurescein- labeled oligonucleotides . 

Cate et a.1. (1991) describe the use of 
oligonucleotide probes directly conjugated to alkaline 
5 phosphatase in combination with a direct chemi luminescent 
substrate (AMPPD) to allow probe detection. 

Labeled probes could readily be purchased form a 
variety of commercial sources, including GENSET, rather then 
synthesized. 

10 Other labels include ligands which can serve as 

specific binding members to a labeled antibody, 
chemiluminescers , enzymes, antibodies which can serve as a 
specific binding pair member for a labeled ligand, and the 
like. A wide variety of labels have been employed in 

15 immunoassays which can readily be employed. Still other 

labels include antigens, groups with specific reactivity, and 
electrochemically detectable moeities. 

In general, labeling of nucleic acids with 
electrophore mass labels ("EML" ) is described, for example, in 

20 Xu et al., J. Chromatography 764.:95-102 (1997). Electrophores 
are compounds that can be detected with high sensitivity by 
electron capture mass spectrometry (EC-MS) . EMLs can be 
attached to a probe using chemistry that is well known in the 
art for reversibly modifying a nucleotide (e.g., well known 

25 nucleotide synthesis chemistry teaches a variety of methods 

for attaching molecules to nucleotides as protecting groups) . 
EMLs are detected using a variety of well known electron 
capture mass spectrometry devices (e.g., devices sold by 
Finnigan Corporation) . Further, techniques that may be used 
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in the detection of EMLs include, for example, fast atomic 
bombardment mass spectrometry (see, e.g., Koster et al . , 
Biomedical Environ. Mass Spec. 14:111-116 (1987)); plasma 
desorption mass spectrometry; electrospray/ionspray (see, 
5 e.g., Fenn et al . , J. Phys. Chem. 88:4451-59 (1984), PCT 

Appln. No. WO 90/14148, Smith et al . , Anal. Chem. 62:882-89 
(1990)); and matrix-assisted laser desorption/ionization 
(Hillenkamp, et al . , "Matrix Assisted UV-Laser 

Desorption/ionization: A New Approach to Mass Spectrometry of 
10 Large Bioraolecules, " Biological Mass Spectrometry (Burlingame 
and McCloskey, eds . ) , Elsevier Science Publishers, Amsterdam, 
pp. 49-60, 1990); Huth-Fehre et al . , "Matrix Assisted Laser 
Desorption Mass Spectrometry of Oligodeoxythymidylic Acids," 
Rapid Communications in Mass Spectrometry, 6:209-13 (1992)) . 

15 In preferred embodmients, the EMLs are attached to a 

probe by a covalent bond that is light sensitive- The EML is 
released from the probe after hybridization with a target 
nucleic acid by a laser or other light source emitting the 
desired wavelength of light. The EML is then fed into a GC-MS 

20 (gas chromatogrph -mass spectrometer) or other appropiate 
device, and identified by its mass. 
EXAMPLE 5 

PrgpaxatilQP Of ggquenqtog Chips an d Arra ys 

A basic example is using 6-mers attached to 50 
25 micron surfaces to give a chip with dimensions of 3 x 3 mm 

which can be combined to give an array of 20 x 20 cm. Another 
example is using 9-mer oligonucleotides attached to 10 x 10 
microns surface to create a 9-mer chip, with dimensions of 5 x 
5 mm. 4000 units of such chips may be used to create a 30 x 
30 30 cm array. In an array in which 4,000 to 16,000 oligochips 
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are arranged into a square array. A plate, or collection of 
tubes, as also depicted, may be packaged with the array as 
part of the sequencing kit. 

The arrays may be separated physically from each 
other or by hydrophobic surfaces. One possible way to utilize 
the hydrophobic strip separation is to use technology such as 
the Iso-Grid Microbiology System produced by QA Laboratories, 
Toronto , Canada . 

Hydrophobic grid membrane filters { HGMF ) have been 
in use in analytical food microbiology for about a decade 
where they exhibit unique attractions of extended numerical 
range and automated counting of colonies. One 

commercially-available grid is ISO-GRID™ from QA Laboratories 
Ltd. (Toronto, Canada) which consists of a square (60 x 60 
cm) of polysulfone polymer (Gelman Tuffryn HT-450, 0.45u pore 
size) on which is printed a black hydrophobic ink grid 
consisting of 1600 (40 x 40) square cells. HGMF have 
previously been inoculated with bacterial suspensions by 
vacuum filtration and incubated on the differential or 
selective media of choice. 

Because the microbial growth is confined to grid 
cells of known position and size on the membrane, the HGMF 
functions more like an MPN apparatus than a conventional plate 
or membrane filter. Peterkin et al . (1987) reported that 
these HGMPs can be used to propagate and store genomic 
libraries when used with a HGMF replicator. One such 
instrument replicates growth from each of the 16 0 0 cells of 
the ISO-GRID and enables many copies of the master HGMF to be 
made (Peterkin et al . , 1987). 
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normalized signals may be scaled, usually from 1-100, to 
compare data from different experiments. Also, in each 
subarray, several control DNAs may be used to determine an 
average background signal in those samples which do not 
5 contain a full match target . For samples obtained from diploid 
(polyploid) scores, homozygotic controls may be used to allow 
recognition of hetero zygotes in the samples. 

H ybridization With Oligonucleotides 

10 Oligonucleotides were either purchased from Genosys 

Inc., Houston, Texas or made on an Applied Biosystems 381A 
DNA synthesizer* Most of the probes used were not purified by 
HPLC or gel electrophoresis. For example, probes were 
designed to have both a single perfectly complementary target 

15 in interferon, a M13 clone containing a 921 bp Eco RI-Bgl II 
human Bl - interferon fragment (Ohno and Tangiuchi, Proc. 
Natl. Acad. flci. 74: 4370-4374 (1981)], and at least one 
target with an end base mismatch in M13 vector itself. 

End labelling of oligonucleotides was performed as 
20 described [Maniatis et al . , Molecular Cloning: A Laboratory 
Manual , Cold Spring Harbor Laboratory Cold Spring Harbor, New 
York (1982)] in 10 ml containing T4 -polynucleotide kinase (5 
units Amersham) , g 32p -ATP (3.3 pM, 10 mCi Amersham 3000 Ci/mM) 
and oligonucleotide (4 pM, 10 ng) . Specific activities of the 
25 probes were 2.5-5 X 10 9 cpm/nM. 

Single stranded DNA (2 to 4 ml in 0.5 NaOH, 1.5 M 
NaCl) was spotted on a Gene Screen membrane wetted with the 
same solution, the filters were neutralized in 0.05 M Na 2 HP0 4 
pH 6,5, baked in an oven at 80°C for 60 min. and UV irradiated 
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emperature. The rate constant £or ^ 

lnCreaseo onl V slightly i„ the o°c to 3 0°c r m » ,„ 

«9«, MoJ . Blo , Mi 361 *J and 

Biol. £2, 383 (197,,, „ „ .„ ' J - 

*- "3(1971,). Hybr ld melting is a first order 
reaction with respect to hybrid concentration (here replaced 

Zr:x to f iiter bouna " ate> - - ■ 

« times t^T e9U " i0n ' *° hytrid "ations 

eb c ancl fc o/ respectively- k i a = -v. 4. 
h , . . n , V iy ' km 18 a ra te constant for 

hybrid melting which is dependent- ™ - 

dependent on temperature and salt 
—ration [Iku ta et a!., M . Acide " 
Por.cli.ce an. Elge n. KoJ . sioJ . ^ " ' 

=1.. J. Mol. Biol. 62; 303 ,1 S71) , n '' &ai9 St 

„h^h ■ . During hybridization 

2 2 a strana aesooiatlon — - -"n g , : r 

strand drsso^atron. reaction takes place ae „eU. Thus the 

amount of hybrid formed i„ time is result of . „ 
-t^^^-h-; result of forward and back 

reaction,. The equilibrium „ay be moved towards hybrid 
for m a tl o» by increasing probe concentration and/or dec! ■ 
temperature. However, during washing cycles! 1 

of buffer, the melting reaction is d ^ 

s reaction ls dominant and the back 
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reaction hybridization is insignificant, since the probe is 
absent. This analysis indicates workable Short Oligonucleotide 
Hybridization (SOH) conditions call be varied for probe 
concentration or temperature. 

5 D or discrimination is defined in equation four: 

D = H p (tj / ^ (t w ) 

H p (t w ) and H ± (t w ) are the amounts hybrids remaining 
after a washing time, t w/ for the identical amounts of 
perfectly and imperfectly complementary duplex, respectively. 
10 For a given temperature, the discrimination D changes with the 
10 length of washing time and reaches the maximal value when H A 
~ B which is equation five. 

The background, B, represents the lowest 
hybridization signal detectable in the system. Since any 

15 further decrease of H ± may not foe examined, D increases upon 
continued washing. Washing past t w just decreases Hp relative 
to B, and is seen as a decrease in D * The optimal washing 
time, t w , for imperfect hybrids, from equation three and 
equation five is: 

20 t w - -In (B / ^ (t 0 ))/ km.i 

Since Hp is being washed for the same t w , combining 
equations, one obtains the optimal discrimination function: 
D = e in (B/Hi < t0)) X H p (t 0 ) / B 

The change of D as a function, of T is important 
25 because of the choice of an optimal washing temperature . It 
is obtained by substituting the Arhenius equation which is: 
K- = Ae ■ Ea ; RT 



into the previous equation to form the final equation: 
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D = H p ( (t 0 )/B X (B/H- (t n (Ap / Ai) e (e 

1 v ° y y a, i 

Wherein B is less than Hi (t 0 ) 



E a,p ,/RT 



Si 



) 



"no, the activation energy £or per£ect 
- ana the activation energy fQr imperfect • „ 

5 either equal, or E leas E D is t 

unctu ii a ^ p d is temperature 
independent, or tersM a, with i nC r ea si„g tel » perature 
respectively. This result ^ ^ ^ 

*i. 9 »t t _ oonditions £or ^ discrlmlM 
son . unjustified. By „ aehing at lower ^^^^ 
obtains egual or better dissipation, nut the time o£ 
washing exponentially increases with che decraaM 
temperature. Discriminate ^ ^ 

if H , . ^imrnat^n more strongly decreases with T, 

ir Hi(t 0 ) increases relative to 1^ (t 0 ) . 

D at lower temperatures depends to a higher degree 
on the H p (t 0 )/ B ratio than on the H (t , / „ ( , f 
m Me v lt , p ^ c o ; / H ± (t 0 ) ratio. 

This result indicates that it is better to oh*, • 
cruantHi-v „<=„•,. better to obtain a sufficient 

quantity of H t m the hybridization regardless of the 
^origination that can be achieved in this step. Better 
d -r^mation can then be obtained by washing, since the 
higher amounts of perfect hybrid allow more time for 
differential melting to snow an effect. Sin , ilarlv , using 
larger amounts of target nudeic acid a necessary 
discrimination can be obtained even with small differences 
between it.,, and K,,,. ° es 

^trapolated to a „ore complex situation than 
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Osi„ g the described theoretical pri 

zt ; ax r riments ' reiiabie J« bl n 

" lth Si * « nudeotides in length M1 

„ere perform „ ith a ^ »1 

Probe, and thus reduced ^ ^ ^ ^ ^ 

hybridizations Th^ hi-„>, 

temperature down to I2°r G ■ • -, 

C ° 12 C - Similarly, the 4-6 X SSC in* „^ ■ 

weuergent m these buffers i« -f^^ ^ ■ . 

^Hty of short oUgonurteotide hybrids was determined 
P-totype octa m er wi th so, GtC contant , J^™* °" * 

TGCTCATG . Th e theoretical expectation is t^^^T"" 

among the less st.hi. . Probe ls 

a xess stable octamers. its t-T-=.r,=.j ^ • 

20 similar to those of m w ^"^zon enthalpy is 

37« <„.,„. Parameter T the! 
so, of the hybrid is _ ltad in Z t ZT I " 

The result shows that T a is ^ c lo „ ! ^ ~ 
25 hybrid th»n f r for tha 8 b P 

hybrid than for an 11 bp duplex [Wallace et ,1 „ , . 

£, 3543 (1 979 , J . al " " UCjelc 

In addition to experiments with model 

oligonucleotides an mi -5 „ A *. 

practical d Ch ° SSn " » "".ten, for , 

30 ! h , d8m ° !l! ' tri,Clon ° f •*« oli goml cleotide 

■5U hybridization. The maim 

6 main Was tO Show uSPfnl • 

w uselul end-mismatch 



15 



0 
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<n™ a ti on „ ith . target siniiar ^ 

™- m „ s a Pplioatlons of tha Mthod Q£ 
way that the M13 vector if qo1f _ . 

base v , ltSelf conta ^s the end mismatched 

base. vector IP , an M13 recombinant contai & 

interferon gene insert ^ • • P n 

9 inSSrt ' Car — ««gle perfectly matched 

r? " 1138 ^^^^ - * ^r nu mber 

of matched targets in comparison fco ^ ^ ^ ^ 

Using low temperature conditions and dot blots, - 
U flC1Snt fences ^ n yb ridi Z ation signals wer e 
between tie dot containing the perfect s „h ^ btaxned 
i . a y Perfect and the mismatched 

; 3 P 7 . an ?" m " — - >— - *> che lar r 

IF-M13 pair of nucleic acids. 

t ar «t hybrldiSaC1 ™ Si ^« 1 *P-d. on the anount of 

-=« Sary oontrol is to 8ho „ that ^^^^^ J - 

. — - the two dota . 8ybridi3ation wtth a - = 

here „ a „ equal amount o£ ^ ^ ^ ^ ^ 
-Kxcx-ncy o£ hybrifl fomatlon increases 
h S13nal £or a dupleJ[ hav±n3 nucieotides 

the mter. Due to their lo „ er moleoular 

~ o f o ligOTUoleo t id e t araet molecul6e =an be ^ 

a t It I: " h * n - — — or „u cl e io 

dCia tnat serves as target. 
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To Masur e the activity of detection with 

ZZT d ?' various ™ of phs - — 

.potted on the fiiter Md hybridl2ed ^ , „ p 

octamer. As little aq rh m-s n ■ 

ie as 50 million unpurifi^H 
5 nr. ™~ ^ "punned phage containina 

method is sufficient p«a«n 4.. 

pr.ctic.lity. ^ ±S Sh ° rt ' - ^ 

As mentioned in the theoretical section above the 
0 equilibrium yield Qf hybrid depends ^ - . he 

I" 7 :;" ° f • - —nee, the sign 

level for the same amount of target with * „ 
. . uar ST e t with 4 nM octamer at l3°r 

is 3 times lower than ™-i*-v, , 

er than w lth a probe concentration of 40 nM and 
decreased 4.5-time* h,, v • ■ ' d 

times by raising the hybridization 
temperature to 25°c. 

The utility of the low temperature wash for 
achieving maximal discrimination -i ^ 

crimination is demonstrated. To make rh* 
Phenomenon visually obvious so M m 

M13 doi- ^ • 50 times more DNA was p ut in the 

M13 dot than m the xp dot using hybridisation with a vector 
specific probe. ln this way , ^ signai ^ 

hybridisation step wit h the actual probe was made strong in 
the, mismatched that in the matched case. The H /H rlT 
1:4. TnvprQ^. * . p /H i ra tio was 

at 7 °c I 7 819,1 in "" Siti - Proles » Mhins 

LhLl : n r d a ratio o£ 2!i - in ™- " - to 

acnieve any discrimination at 25 V «^ 

si£Jnal . n at 25 c, since the matched target 

aaaxng, at the same tnm*=> hi^~ 

cime, the signal from the* 
hatched h ybrid ia still ^ ^ 

^c^inatio„ at ^ COWed to ^ ^ m so bM ^ 
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clearly visi bl e. „ one oonBlders the so mimte po±nt ^ ^ 

and the « minut9 point , t 13 „ c ^ ^ 

-gnal near the laval , „ hich represents 

times for the respective condicions> ifc is ob J us 

^ cnat the amount of qp^orai - ■ 

u oi: several times greater at 7°r «-u 0 

yx«dter at 7 C than at 13°C. 

To illustrate this further, the tlM bourse of ^ 
^crimination „ lth WMhing o£ same Mount ^ 

Vhr.a at the two temperatures shows the higher maxirai » 

he lower temperature. ThMe Msulte confirm ^ 
hange of o with temperature ^ the ratlo of Mounts ■» 
t»o types of hybrld at the etarc of w5Bh . n3 ^ 

olico , In ° rder " Sh °" 9Cneral UtUity ° £ tte *»~t 
^ * conditions. we have Io o k ed 

^ hyhrrd^ation o f 4 heptamers. „ octaMrs ^ m 

probes up to „ nucieotiaes in l ength in our simple M13 

Z2Zr Th " e in ° 1Ude - tte — « -a oct.mer 

GGCaGGCG resenting the two extremes of ac COI5 tent 

it::* sc f ™ and sequenee ~ ^* - «- 

stahrl.ty of short hyhrias IBMSalMet « al ., ^ 

short oligonucleotide conditions were ap pl i 0 a«,l. to all testea 
P~h.es » achieving su£fioient aisorimlnation . 
^cr^nation value obCalMd „ lth probes ^ nu=ieo 

length was 20, a several fvvi/j * 
5 ±€ , f ° ld dr ° P due to sequence variation 

5 is easily tolerated. 

The M!3 syst™ has the advantage of showing the 
effects of target DBA complexity on the levels of 
discrimination. F or two octamers having either none or five 
■^matched targets ,„a aiffering in only one GC pair the 
ohservea discriminations were ia.3 ana respectively 
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In order to show tha utility of this method, three 
probes 8 nucleotides i„ l e „ g th „ ere tested on q£ 

plasm.d DNA dots made from . library Blue3<=ript ve 

one probe was present and specific tor adscript vector but ' 
was absent in M!3, while the other two probes had targets that 
were inserts o f knOTm seance. This ^ ^ 
Of hybridisation negative or positive controi DNAs with each 
probe. This probe sequence (CTCCCTTT) also had a 
complementary target in th . interferon insert. since th e M13 
dot „ negatrve „ hile the interferon insert in either «3 or 
Bluescr.pt „ aa positive, the hybridisation is sequence 
specific. Similarly, probes that detect th. tar g et seance 
n only of 51 insartei or ^ none ^ exsmined 

alo„ g „ lt h controls that confirm that hybridization would have 
15 occurred ,f the appropriate targets were present in the 

clones. 

Thermal stability curves for very short 

aTilT^T " e that are 6 - 8 in — - 

t leas is c lower than for hybrids al - 12 nucleotides in 
0 length [Flg . , a nd Wallace et al . . Meio Aalds Res 

3S.3-35S, ( „ 7S)J . Sowever, performing the hybridation 
r Ctl ° n " S - -ith > very practical 

nM concentration of ol ig onucleotiae probe allows the detection 
of complementary seance in a kno „„ or ^ ^ ^ 
target. To determine an unknown nucleic acid seguence 
completely, an entire set containing 65 ,S3 5 s-mer probes may 
- used, sufficient amounts of nucleic acid for this purpose 
a e present in convenient biological samples such as a \Z 
m.crolrters of M 13 culture, a piasmid prep from iO ml of 
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bacterial culture or a single colony of bacteria or l 

, i Mtceria, or less than 

1 ml of a standard PCR reaction. 

Short oligonucleotides 6 -l 0 nucleotides long give 
excellent discrimination. The relate * 
5 ^k-t-- • x. ^latxve decrease in hybrid 

stability W1 th a single end mismatch is great£r than ^ 

longer probes. Kesults with the octamer TOCTCA.G support this 
conclusion. In the experiments, the target with a G/T end 
mismatch, hybridization to the target of this type of 
mismatch is the most stable of all other types of 
' oligonucleotide. This discrimination achieved is the same as 

J r 6 "^ than int6rnal ^ * * » *~ paire d 

UP 6X » O/T mismatch in a „ P L ed 

ex [Ikuta et al ., Nucl> Acids _ P 

Exploiting these discrimination properties using the described 
hybridation conditions for short oligonucleotide 
hybridization allows a very precise determination of 
oligonucleotide targets. In contrast to the ease of detecting 

discrimination between perfect ^ • * 9 

Perrect and imperfect hybrids a 
problem that mav exi^t- • ' 

mcty exist with using verv ^hor-t- ~n ■ 
. y ver y snort oligonucleotides 

x- the preparation of sufficient amounts of hybrids In 
P-ctice, the need to discriminate Hp and Hi is aided by 
increasing the amount of DNA in the dot and/or the probe 

^ t dCCreaSing ™-ion -mperature. 

However, higher probe concentrations usually increase 

background. Moreover, there are limits to ^ 

«xe limits to the amounts of 
target nucleic acid that are m-art-^i - 

die practical to use T>n - « v>-v-~i~t 
~ solved by the hlgher oonoantration of thea J> - I — 

Probe. Pu rther improvements ray be e 

of ~„pe titoro £or unspeoi£i= binding ^ ^ ^-J-« 
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by changing the hybridization support material. Moreover, for 
probes having E a less than 45 Kcal/mol (e.g. for many heptamers 
and a majority of hexamers, modified oligonucleotides give a 
more stable hybrid [Asseline, et al . , Proa. Nat'l Acad. Sci. 
81 : 3297 (1984)] than their unmodified counterparts. The 
hybridization conditions described in this invention for short 
oligonucleotide hybridization using low temperatures give 
better discriminating for all sequences and duplex hybrid 
inputs. The only price paid in achieving uniformity in 
hybridization conditions for different sequences is an 
increase in washing time from minutes to up to 24 hours 
depending on the sequence. Moreover; the washing time can be 
further reduced by decreasing the salt concentration. 

Although there is excellent discrimination of one 
matched hybrid over a mismatched hybrids, in short 
oligonucleotide hybridization, signals from mismatched hybrids 
exist, with the majority of the mismatch hybrids resulting 
from end mismatch. This may limit insert sizes that may be 
effectively examined by a probe of a certain length. 

The influence of sequence complexity on 
discrimination cannot be ignored. However, the complexity 
effects are more significant when defining sequence 
information by short oligonucleotide hybridization for 
specific, nonrandom sequences, and can be overcome by using an 
appropriate probe to target length ratio. The length ratio is 
chosen to make unlikely, on statistical grounds, the 
occurrence of specific sequences which have a number of 
end-mismatches which would be able to eliminate or falsely 
invert discrimination. Results suggest the use of 
oligonucleotides 6, 7, and 8 nucleotides in length on target 
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nucleic acid inserts shorter than 0.6, 2.5, and 10 kb, 
respectively. 

EXAMPT.E! I-; 
DNA St^qi^r.^ 

An array of subarrays allows for efficient 
sequencing of a small set of samples arrayed in the form of 
replicated subarrays; Per example, 64 samples may bg arrayed 
on a a x 8 mm subarray and 16 x 24 sub arrays may be replicated 
on a 15 X 23 cm membrane with l mm wide spacers between ^ 
subarrays. Several replica membranes may be made. For example 
probes from a universal set of three thousand seventy-two ' 
7-mers may be divided in thirty-two 96-well plates and 
labelled by kinasing. Pour membranes may be processed in 
Parallel during one hybridation cycle. On each membrane, 384 
probes may be scored. All probes may be scored in two 
hybridization cycles. Hybridisation intensities may be scored 
and the sequence assembled as described below. 

If a single sample subarray or subarrays contains 
several unknowns, especially when similar samples are used, a 
smaller number of probes may be sufficient if the y are 
mtelligently selected on the basis of results of previously 
scored probes. Por example, if probe ±s ^ fc 

there is a small chance that any of 8 ov.n^' 

or o overlapping probes are 

pos.tive. Ths sconcing process in this case oonBists Qf 

7' • <- -inl^lly ove rl ,p pad probes to 

aetme positive anchors anr3 i-h^r, <- 

and then to successively select probes 
which confirms one of the mocst- hi, -\ v 

tne most likely hypotheses about the 
order of anchors and size and t-™^ * 

and ty P e of ^aps between them. m 
this second phase, nool=s o-f o -i n 

e, pools of 2-10 probes may be used where each 
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probe is selected to be pos±c . ve , n oniy Dm sampie ^ 

different from the sables expected to be positive with 

other probes from the pool. 

The subarray approach allows efficient 
5 implementation of probe competition (overlapped probes) or 
probe cooperation (continuous stacking of probes) in solving 
branching probiems . After hybridation of a universal set of 
probes the sequence assembly program determines candidate 
sequence subfragments ( SFs ) . Por the further assembly of SPs 
additional information has to be provided (from overlapped ' 

sequences of DNA fragments, si mi i ar seque nces Bind 

ocyuences, single pass crel 

sequences, or from other hybridization or restriction m ■ 
. H=1+ . „ . resurxction mapping 

data) . Primers f or single pass seguenc . ^ fchrough fche 

branch points are identified from the* c-mr «, 

w - m the SBH 8e quence information 

or from known veccor eeguences _ e 3 _ fiankins 

to the vector insert site. and standard Sanger-sequencing 
reactions are perform on the sample DNA. The sequence 
obtained from this single pass gel sequencing is compared to 
the Sfs that read into and out of the branch points to 
identify the order of the sfs. Purther, singel pass gol 

sequencing may be combined with SBB i-„ h. 

"° iBH to de novo sequence or re- 
sequence a nucleic acid. 

Competitive hybridization and continuous stacking 
interactions can also be used to assemMe S£s . TheS e 
approaches are of Umited commercial value for sequencing of 
large numbers of samples by SBH wherein a labelled probe is 
applied to a sample affixed to an array if . „ nl£or m array is 

rep lie, r Mtely ' * MlySi3 ° f 3ma11 nUmte " ° f ^ 
repuca subarrays allows efficient implementation of both 

approaches, on each cf the replica subarrays, one branching 
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point may be tested for one or more DNA samples using pools of 
probes similarly as in solving mutated sequences in different 
samples spotted in the same subarray (see above) . 

If in each of 64 samples described in this example, 
5 there are about 100 branching point s, and if 8 samples are 
analyzed in parallel in each subarray, then at least 800 
subarray probings solve all branches. This means that for the 
3072 basic probings an additional 800 probings (25%) are 
employed. More preferably/ two probings are used for one 

10 branching point. If the subarrays are smaller, less 

additional probings are used. For example, if subarrays 
consist of 16 samples, 200 additional probings may be scored 
(6%) . By using 7-mer probes (N 1 , 2 B 7 N 1 _ 2 ) and competitive or 
collaborative branching solving approaches or both, fragments 

15 of about 1000 bp fragments may be assembled by about 4 00 0 
probings. Furthermore, using 8-mer probes (NB a N) 4 kb or 
longer fragments may be assembled with 12,000 probings. 
Gapped probes, for example, NB 4 NB 3 N or NB 4 NB 4 N may be used to 
reduce the number of branching points. 

20 BXMfl P Mi 12 

PNA Analysis bv Transient Attachment to Subarrays 
of Probes and Ligation of Labelled Probes 



Oligonucleotide probes having an informative length 
25 of four to 40 bases are synthesized by standard chemistry and 
stored in tubes or in multiwell plates. Specific sets of 
probes comprising one to 10,000 probes are arrayed by 
deposition or in situ synthesis on separate supports or 
distinct sections of a larger support. In the last case, 
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sections or subarrays may be sep arated by physical or 
hydrophobic barriers. The probe arrays may be prepared by ia 
situ synthesis. A sample DNA of appropriate size is hybridized 
with one or more specific arrays. Many samples may be 
interrogated as pools at the same subarrays or independently 
with different subarrays within one support. Simultaneously 
with the sample or subsequently, a single labelled probe or a 
Pool of labelled probes is added on each of the subarrays if 
attached and labelled probes hybridize back to back on the 
complementary target in the sample DNA fchey ^ 

Occurrence of ligation will be me asured by detecting a label 
from the probe. 

This procedure is a variant of the descr ibed DNA 
analysis process in which DNA samples are not permanently 
attached to the support. Transient attachment is provided by 
Probes fixed to the support, m this case there i s no need for 
a target DNA arraying process. In addition, Ugation allows 
detection of longer oligonucleotide fences by combining 
short labelled probes with short fixed probes. 

The process has several unique features. Basically 
the transient attachment of the target allows its reuse After 
ligation occur the target may be released and the label will 
stay covalently attached to the support. This feature allows 
cycling the target and production of detectable signal with a 
-all quantity of the target. Under optimal conditions, 
targets do not need to be ampli fied , e . g . nafcural ^ 
the DNA samples may be directly used for diagnostics and 
sequencing purposes. Targets may be released by cycling the 
temperature between efficient hybridization and efficient 
meltmg of duplexes. More preferablly, there is no cycling 
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The temperature and concentrations of components may be 
defined to have an equilibrium between free targets and 
targets entered in hybrids at about 50:50% level. In this case 
there is a continuous production of ligated products. For 
5 different purposes different equilibrium ratios are optimal . 

An electric field may be used to enhance target use. 
At the beginning; a horizontal field pulsing within each 
subarray may be employed to provide for faster target sorting. 
In this phase, the equilibrium is moved toward hybrid 

10 formation, and unlabelled probes may be used. After a target 
sorting phase, an appropriate washing (which may be helped by 
a vertical electric field for restricting movement of the 
samples) may be performed. Several cycles of discriminative 
hybrid melting, target harvesting by hybridization and 

15 ligation and removing of unused targets may be introduced to 
increase specificity. In the next step, labelled probes are 
added and vertical electrical pulses may be applied. By 
increasing temperature, an optimal free and hybridized target 
ratio may be achieved. The vertical electric field prevents 

20 diffusion of the sorted targets. 

The subarrays of fixed probes and sets of labelled 
probes (specially designed or selected from a universal probe 
set) may be arranged in various ways to allow an efficient and 
flexible sequencing and diagnostics process. For example, if a 

25 short fragment (about 100-500 bp) of a bacterial genome is to 
be partially or completely sequenced, small arrays of probes 
(5-30 bases in length) designed on the bases of known sequence 
may be used. If interrogated with a different pool of 10 
labelled probes per subarray, an array of 10 subarrays each 

30 having 10 probes, allows checking of 200 bases, assuming that 
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only two bases connected by ligation are scored. Under the 
conditions where mismatches are discriminated throughout the 
hybrid, probes may be displaced by more than one base to cover 
the longer target with the same number of probes. By using 
5 long probes, the target may be interrogated directly without 
amplification or isolation from the rest of DNA. in the sample. 
Also, several targets may be analyzed (screened for) in one 
sample simultaneously. If the obtained results indicate 
occurrence of a mutation (or a pathogen) , additional pools of 
10 probes may be used to detect type of the mutation or subtype 
of pathogen. This is a desirable feature of the process which 
may be very cost effective in preventive diagnosis where only 
a small fraction of patients is expected to have an infection 
or mutation. 

15 In the processes described in the examples, various 

detection methods may be used, for example, radiolabels, 
fluorescent labels, enzymes or antibodies (chemi luminescence) , 
large molecules or particles detectable by light scattering or 
interf erometric procedures* 

20 E XAM P LE 3-3 

Secruenci^cr a Target XTsing Qctamers and Nonamers 

Data resulting from the hybridization of octamer and 
nonamer oligonucleotides shows that sequencing by 
hybridizat ion provides an extremely high degree of accuracy. 
25 In this experiment, a known sequence was used to predict a 

series of contiguous overlapping component octamer and nonamer 
oligonucleotides . 



In addition to the perfectly matching 
oligonucleotides, mismatch oligonucleotides, mismatch 
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oligonucleotides wherein internal or end mismatches occur in 
the duplex formed by the oligonucleotide and the target were 
examined. In these analyses, the lowest practical temperature 
was used to maximize hybridization formation. Washes were 
5 accomplished at the same or lower temperatures to ensure 

maximal discrimination by utilizing the greater dissociation 
rate of mismatch versus matched oligonucleotide/target 
hybridization. These conditions are shown to be applicable to 
all sequences although the absolute hybridization yield is 
10 shown to be sequence dependent . 

The least destabilizing mismatch that can be 
postulated is a simple end mismatch, so that the test of 
sequencing by hybridization is the ability to discriminate 
perfectly matched oligonucleotide/target duplexes from 
15 end -mismatched oligonucleotide/target duplexes. 

The discriminative values for 102 of 105 hybridizing 
oligonucleotides in a dot blot format were greater than 2 
allowing a highly accurate generation of the sequence. This 
system also allowed an analysis of the effect of sequence on 
20 hybridization formation and hybridization instability. 

One hundred base pairs of a known portion of a 
human- interferon genes prepared by PCR # i.e. a 10 0 bp target 
sequence, was generated with data resulting from the 
hybridization of 105 oligonucleotides probes of known sequence 
25 to the target nucleic acid. The oligonucleotide probes used 
included 72 octamer and 21 nonatner oligonucleotides whose 
sequence was perfectly complementary to the target. The set 
of 93 probes provided consecutive overlapping frames of the 
target sequence e displaced by one or two bases . 
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To evaluate the effect of mismatches, hybridization 
was examined for 12 additional probes that contained at least 
one end mismatch when hybridized to the 100 bp test target 
sequence. Also tested was the hybridization of twelve probes. 
5 with target end-mismatched to four other control nucleic acid 
sequences chosen so that the 12 oligonucleotides formed 
perfectly matched duplex hybrids with the four control DNAs . 
Thus, the hybridization of internal mismatched, end-mismatched 
and perfectly matched duplex pairs of oligonucleotide and 

10 target were evaluated for each oligonucleotide used in the 

experiment. The effect of absolute DNA target concentration on 
the hybridization with the test octamer and nonamer 
oligonucleotides was determined by defining target DNA 
concentration by detecting hybridization of a different 

15 oligonucleotide probe to a single occurrence non- target site 
within the co-amplified plasmid DNA. 

The results of this experiment showed that all 
oligonucleotides containing perfect matching complementary 
sequence to the target or control DNA hybridized more strongly 

20 than those oligonucleotides having mismatches. To come to 
this conclusion, we examined H p and D values for each probe. 
H p defines the amount of hybrid duplex formed between a test 
target and an oligonucleotide probe. By assigning values of 
between 0 and 10 to the hybridization obtained for the 105 

25 probes, it was apparent that 68.5% of the 105 probes had an H p 
greater than 2 . 

Discrimination (D) values were obtained where D was 

defined as the ratio of signal intensities between 1) the dot 
containing a perfect matched duplex formed between test 

30 oligonucleotide and target or control nucleic acid and 2) the 
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dot containing a mismatch duplex formed between the same 
oligonucleotide and a different site within the target or 
control nucleic acid. Variations in the value of D result 
from either 1) perturbations in the hybridization efficiency 
5 which allows visualization of signal over background, or' 2) 
the type of mismatch found between the test oligonucleotide 
and the target. The D values obtained in this experiment were 
between 2 and 40 for 102 of the 105 oligonucleotide probes 
examined. Calculations of D for the group of 102 
10 oligonucleotides as a whole showed the average D was 10.6. 

There were 20 cases where oligonucleotide/target 
duplexes exhibited an end-mismatch. In five of these, D was 
greater than 10. The large D value in these cases is most 
likely due to hybridization destabilization caused by other 
15 than the most stable (G/T and G/A) end mismatches. The other 
possibility is there was an error in the sequence of either 
the oligonucleotides or the target. 

Error in the target for probes with low H p was 
excluded as a possibility because such an error would have 

20 affected the hybridization of each of- the other eight 
overlapping oligonucleotides. There was no apparent 
instability due to sequence mismatch for the other overlapping 
oligonucleotides, indicating the target sequence was correct. 
Error in the oligonucleotide sequence was excluded as a 

25 possibility after the hybridization of seven newly synthesized 
oligonucleotides was re-examined. Only 1 of the seven 
oligonucleotides resulted in a better D value. Low hybrid 
formation values may result from hybrid instability or from an 
inability to form hybrid duplex. An inability to form hybrid 

30 duplexes would result from either 1) self complementarity of 
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the chosen probe or 2) target/target self hybridization. 
Oligonucleotide/oligonucleotide duplex formation may be 
favored over oligonucleotide/target hybrid duplex formation if 
the probe was self -complementary . Similarly, target/target 
5 association may be favored if the target was 

self -complementary or may form internal palindromes. In 
evaluating these possibilities, it was apparent from probe 
analysis that the questionable probes did not form hybrids 
with themselves. Moreover, in examining the contribution of 

10 target/target hybridization, it was determined that one of the 
questionable oligonucleotide probes hybridized inefficiently 
with two different DNAs containing the same target . The low 
probability that two different DNAs have a self -complementary 
region for the same target sequence leads to the conclusion 

15 that target /target hybridization did not contribute to low 
hybridization formation. Thus, these results indicate that 
hybrid instability and not the inability to form hybrids was 
the cause of the low hybrid formation observed for specific 
oligonucleotides. The results also indicate that low hybrid 

20 formation is due to the specific sequences of certain 
oligonucleotides. Moreover, the results indicate that 
reliable results may be obtained to generate sequences if 
octamer and nonamer oligonucleotides are used. 

These results show that using the methods described 
25 long sequences of any specific target nucleic acid may be 
generated by maximal and unique overlap of constituent 
oligonucleotides. Such sequencing methods are dependent on 
the content of the individual component oligomers regardless 
of their frequency and their position. 
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The sequence which is generated using the algorithm 
described below is of high fidelity. The algorithm tolerates 
false positive signals from the hybridization dots as is 
indicated from the fact the sequence generated from the 105 
5 hybridization values, which included four less reliable 
values, was correct. This fidelity in sequencing by 
hybridization is due to the "all or none" kinetics of short 
oligonucleotide hybridization and the difference in duplex 
stability that .exists between perfectly matched duplexes and 

10 mismatched duplexes. The ratio of duplex stability of matched 
and end-mismatched duplexes increases with decreasing duplex 
length. Moreover, binding energy decreases with decreasing 
duplex length resulting in a lower hybridization efficiency. 
However, the results provided show that octamer hybridization 

15 allows the balancing of the factors affecting duplex stability 
and discrimination to produce a highly accurate method of 
sequencing by hybridization. Results presented in other 
examples show that oligonucleotides that are 6, 7, or 8 
nucleotides can be effectively used to generate reliable 

20 sequence on targets that are 0 . 5 kb (for hexamers) 2 kb (for 
septamers) and 6kb (for octamers) . The sequence of long 
fragments may be overlapped to generate a complete genome 
sequence , 
ex am p le 14 

25 Analyzing the Data Obta i n e d 

Image files are analyzed by an image analysis 
program, like DOTS program (Drmanac et al . , 1993), and scaled 
and evaluated by statistical functions included, e.g., in 
SCORES program (Drmanac et al . 1994). From the distribution 
30 of the signals an optimal threshold is determined for 
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transforming signal into +/- output. From the position of the 
label detected, F + P nucleotide sequences from the fragments 
would be determined by combining the known sequences of the 
immobilized and labeled probes corresponding to the labeled 
5 positions* The complete nucleic acid sequence or sequence 
subf ragments of the original molecule, such as a human 
chromosome, would then be assembled from the overlapping F + P 
sequence determined by computational deduction. 

One option is to transform hybridization signals 

10 e.g., scores, into +/- output during the sequence assembly 
process. In this case, assembly will start with a F +■ P 
sequence with a very high score, for example F + P sequence 
AAAAAATTTTTT . Scores of all four possible overlapping probes 
AAAAATTTTTTA , AAAAATTTTTT , AAAAATTTTTTC and AAAAATTTTTTG 

15 and three additional probes that are different at the 

beginning ( TAAAAATTTTTT , ; CAAAAATTTTTT , ; GAAAAATTTTTT, are 
compared and three outcomes defined: (i) only the starting 
probe and only one of the four overlapping proves have scores 
that are significantly positive relatively to the other, six 

20 probes, in this case the AAAAAATTTTTT sequence will be 

extended for one nucleotide to the right; (ii) no one probe 
except the starting probe has a significantly positive score, 
assembly will stop, e.g., the AAAAAATTTTTT sequence is at the 
end of the DNA molecule that is sequenced; (iii) more than one 

25 significantly positive probe among the overlapped and/or other 
three probes is found; assembly is stopped because of the 
error or branching (Drmanac et al . , 1989). 

The processes of computational deduction would 
employ computer programs using existing algorithms (see, 
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e.g., Pevzner, 1989/ Drmanac et al . , 1991; Labat and Drmanac, 
1993; each incorporated herein by reference). 

If, in addition to F + P, F (space 1)P, F (space 
2)P, F( space 3)P or F (space 4)P are determined, algorithms 
will be used to match all data sets to correct potential 
errors or to solve the situation where there is a branching 
problem (see, e.g., Drmanac et al - , 1989; Bains et al. , 1988; 
each incorporated herein by reference) . 

Conducting Sequencing bv Two Step Hybridization 

Following the certain examples to describe the 
execution of the sequencing methodology contemplated by the 
inventor. First, the whole chip would be hybridized with 
mixture of DNA as complex as 100 million of bp (one human 
chromosome) . Guidelines for conducting hybridization can be 
found in papers such as Drmanac et al. (1990); Khrapko et al . 
(1991); and Broude et al . (1994) . These articles teach the 
ranges of hybridization temperatures, buffers and washing( 
steps that are appropriate for use in the initial steps of 
Format 3 SBH. 

The present inventor particularly contemplates that 
hybridization is to be carried out for up to several hours in 
high salt concentrations at a low temperature (-2°C to 5°C) 
because of a relatively low concentration of target DNA that 
can be provided. For this purpose, SBC buffer is used instead 
of sodium phosphate buffer (Drmanac et al , 1990) , which 
precipitates at 10°C. Washing does not have to be extensive (a 
few minutes) because of the second step; and can be completely 
eliminated when the hybridization cycling is used for the 
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sequencing of highly complex DNA samples. The same buffer is 
used for hybridization and washing steps to be able to 
continue with the second hybridization step with labeled 
probes . 

5 After proper washing using a simple robotic device 

on each array, e.g., a 8 x 8 mm array, one labeled, probe, 
e.g., a 6-mer, would be added. A 96-tip or 96-pin device 
would be used, performing this in 42 operations. Again, a 
range of discriminatory conditions could be employed, as 
10 previously described in the scientific literature. 

The present inventor particularly contemplates the 
use of the following conditions. First, after adding labeled 
probes and incubating for several minutes only (because of the 
high concentration of added oligonucleotides) at a low 

15 temperature (0-5°C) , the temperature is increased to 3-lO°C, 
depending on F + P length, and the washing buffer is added. 
At this time, the washing buffer used is one compatible with 
any ligation reaction (e.g., 100 raM salt concentration range). 
After adding ligase, the temperate is increased again to 

20 15~37°C to allow fast ligation (less than 30 min) and further 
discrimination of full match and mismatch hybrids. 

The use of cat ionic detergents is also contemplated 
for use in Format 3 SBH, as described by Pontius & Berg (19 91, 
incorporated herein by reference) . These authors describe the 
25 use of two simple cat ionic detergents, dodecy- and 

cetyltrimethylammonium bromide (DTAB and CTAB) in DNA 
renaturation . 



DTAB and CTAB are variants of the quaternary amine 
tetramethylammonium bromide (TMAB) in which one of the methyl 
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groups is replaced by either a 12 -carbon (DTAB) or a 16-carbon 
(CTAB) alkyl group, TMAB is the bromide salt of the 
tetramethylamrnonium ion, a reagent used in nucleic acid 
renaturation experiments to decrease the G-C content bias of 
5 the melting temperature, DTAB and CTAB are similar in 
structure to sodium dodecyl sulfate (SDS) , with the 
replacement of the negatively charged sulfate of SDS by a 
positively charged quaternary amine. While SDS is commonly 
used in hybridization buffers to reduce nonspecific binding 
10 and inhibit nucleases, it does not greatly affect the rate of 
renaturation . 

When using a ligation process, the enzyme could be 
added with the labeled probes or after the proper washing step 
to reduce the background. Although not previously proposed 

15 for use in any SBH method, ligase technology is well 

established within the field of molecular biology. For 
example, Hood and colleagues described a ligase -mediated gene 
detection technique (Landegren et al . , 1988), the methodology 
of which can be readily adapted for use in Format 3 SBH. Wu 

20 & Wallace also describe the use of bacteriophage T4 DNA ligase 
to join two adjacent, short synthetic olignucleotides . Their 
oligo ligation reactions were carried out in 50 mM Tris HCl pH 
7.6, 10 mM MgCl 2 , 1 mM ATP, 1 mM DTT, and 5% PEG. Ligation 
reactions were heated to 10 0 P C for 5-10 min followed by cooling 

25 to 0°C prior to the addition of T4 DNA ligase (1 unit/ Bethesda 
Research Laboratory) . Most ligation reactions were carried 
out at 30°C and terminated by heating to 100°C for 5 min. 



30 



Final washing appropriate for discriminating 
detection of hybridized adjacent, or ligated, oligonucleotides 
of length (F + P) , is then performed. This washing step is 
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done in water for several minutes at 4 0-60°C to wash out all 
the non-ligated labeled probes, and all other compounds, to 
maximally reduce background. Because of the covalently bound 
labeled oligonucleotides , detection is simplified (it does not 
5 have time and low temperature constrains) . 

Depending on the label used, imaging of the chips is 
done with different apparati. For radioactive labels, 
phosphor storage screen technology and Phosphorlmager as a 
scanner may be used (Molecular Dynamics, Sunnyvale, CA) . 

10 Chips are put in a cassette and covered by a phosphorous 

screen. After 1-4 hours of exposure, the screen is scanned 
and the image file stored at a computer hard disc. For the 
detection of fluorescent labels, CCD cameras and 
epif luorescent or confocal microscopy are used. For the chips 

15 generated directly on the pixels of a CCD camera, detection 
can be performed as described by Eggers et al . (1994, 
incorporated herein by reference) . 

Charge -coup led device (CCD) detectors serve as 
active solid supports that quantitatively detect and image the 

20 distribution of labeled target molecules in probe-based 

assays. These devices use the inherent characteristics of 
microelectronics that accommodate highly parallel assays, 
ultrasensitive detection, high throughput, integrated data 
acquisition and computation. Eggers et al . (1994) describe 

25 CCDs for use with probe-based assays, such as Format 3 SBH of 
the present invention, that allow quantitative assessment 
within seconds due to the high sensitivity and direct coupling 
employed. 
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The integrated CCD detection approach enables the 
detection of molecular binding events on chips. The detector 
rapidly generates a two-dimensional pattern that uniquely 
characterizes the sample. In the specific operation of the 
5 CCD-based molecular detector, distinct biological probes are 
immobilized directly on the pixels of a CCD or can be attached 
to a disposable cover slip placed on the CCD surface. The 
sample molecules can be labeled with radioisotope, 
chemiluminescent or fluorescent tags, 

10 Upon exposure of the sample to the CCD-based probe 

array, photons or radioisotope decay products are emitted at 
the pixel locations where the sample has bound, in the case of 
Format 3, to two complementary probes. In turn, electron-hole 
pairs are generated in the silicon when the charged particles, 

15 or radiation from the labeled sample, are incident on the CCD 
gates. Electrons are then collected beneath adjacent CCD 
gates and sequentially read out on a display module. The 
number of photoelectrons generated at each pixel is directly 
proportional to the number of molecular binding events in such 

20 proximity. Consequently, molecular binding can be 
quantitatively determined (Eggers et al . , 1994), 

By placing the imaging array, in proximity to the 
sample, the collection efficiency is improved by a factor of 
at least 10 over lens-based techniques such as those found in 
25 conventional CCD cameras. That is, the sample (emitter) is in 
near contact with the detector (imaging array) , and this 
eliminates conventional imaging optics such as lenses and 
mirrors . 
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When radioisotopes are attached as reporter groups 
to the target molecules, energetic particles are detected. 
Several reporter groups that emit particles of varying 
energies have been successfully utilized with the 
5 micro -fabricated detectors, including 32 P, 33 P, 35 S, 14 C and 12S L. 
The higher energy particles, such as from 32 P / provide the 
highest molecular detection sensitivity, whereas the lower 
energy particles, such as from 35 S, provide better resolution. 
Hence the choice of the radioisotope reporter can be tailored 
10 as required. Once the particular radioisotope label is 
selected, the detection performance can be predicted by 
calculating the signal-to-noise ration (SNR) , as described by 
Sggers et al . (1994). 

An alternative luminescent detection procedure 
15 involves the use of fluorescent or chemi luminescent reporter 
groups attached to the target molecules. The fluorescent 
labels can be attached covalently or through interaction. 
Fluorescent dyes, such as ethidium bromide, with intense 
absorption bands in the near UV (300-350 nm) range and 
20 principal emission bands in the visible (500-650 nm) range, 

are most suited for the CCD devices employed since the quantum 
efficiency is several orders of magnitude lower at the 
excitation wavelength then at the fluorescent signal 
wavelength . 

25 Prom the perspective of detecting luminescence, the 

polysilicon CCD gates have the built-in capacity to filter 
away the contribution of incident light in the UV range, yet 
are very sensitive to the visible luminescence generated by 
the fluorescent reporter groups. Such inherently large 

30 discrimination against UV excitation enables, large SNRs 
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(greater than 10 0) to be achieved by the CCDs as formulated in 
the incorporated paper by Eggers et al . (1994) . 

For probe immobilization on the detector, 
hybridization matrices may be produced on inexpensive Si0 2 
5 wafers, which are subsequently placed on the surface of the 
CCD following hybridization and drying. This format is 
economically efficient since the hybridization of the DNA is 
conducted on inexpensive disposable Si0 2 wafers, thus allowing 
reuse of the more expensive CCD detector. Alternatively, the 
10 probes can be immobilized directly on the CCD to create a 
dedicated probe matrix. 

To immobilize probes upon the Si0 2 coating, a uniform 
epoxide layer is linked to the film surface, employing an 
epoxy-silane reagent and standard Si0 2 modification chemistry. 

15 Amine-modif ied oligonucleotide probes are then linked to the 
Si0 2 surface by means of secondary amine formation with the 
epoxide ring. The resulting linkage provides 17 rotatable 
bonds of separation between the 3 base of the oligonucleotide 
and the Si0 2 surface. To ensure complete amine deprotonation 

20 and to minimize secondary structure formation during coupling, 
the reaction is performed in 0.1 M KOH and incubated at 37°C 
for 6 hours. 

In Format 3 SBH in general, signals are scored per 
each of billion points. It would not be necessary to 
25 hybridize all arrays, e.g., 4000 5x5 mm, at a time and the 
successive use of smaller number of arrays is possible. 

Cycling hybridizations are one possible method for 
increasing the hybridization signal. In one cycle, most of 
the fixed probes will hybridize with DNA fragments with tail 
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sequences non-complementary for labeled probes. By increasing 
the temperature, those hybrids will be melted. In the next 
cycle , some of them (-0.1%) will hybridize with an 
appropriate DNA fragment and additional labeled probes will be 
5 ligated. In this case, there occurs a discriminative melting 
of DNA hybrids with mismatches for both probe sets 
simultaneously. 

In the cycle hybridization, all components are added 
before the cycling starts, at the 3 7°C for T4 , or a higher 

10 temperature for a thermostable ligase. Then the temperature 
is decreased to 15-37°C and the chip is incubated for up to 10 
minutes, and then the temperature is increased to 3 7°C or 
higher for a few minutes and then again reduced. Cycles can 
be repeated up to 10 times. In one variant, an optimal higher 

15 temperature (10-50 P C) can be used without cycling and longer 
ligation reaction can be performed (1-3 hours) . 

The procedure described herein allows complex chip 
manufacturing using standard synthesis and precise spotting of 
oligonucleotides because a relatively small number of 
20 oligonucleotides are necessary. For example, if all 7~mer 
oligos are synthesized (16384 probes) , lists of 256 million 
14-mers can be determined. 

One important variant of the invented method is to 
use more than one differently labeled probe per base array. 
25 This can be executed with two purposes in mind; multiplexing 
to reduce number of separately hybridized arrays; or to 
determine a list of even longer oligosequences such as 3 x 6 
or 3 x 7 . In this case, if two labels are used, the 
specificity of the 3 consecutive oligonucleotides can be 
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almost absolute because positive sites must have enough 
signals of both labels. 

A further and additional variant is to use chips 
containing BxNy probes with y being from 1 to 4 . Those chips 
5 allow sequence reading in different frames. This can also be 
achieved by using appropriate sets of labeled probes or both F 
and P probes could have some unspecified end positions (i.e., 
some element of terminal degeneracy) . Universal bases may 
also be employed as part of a linker to join the probes of 
10 defined sequence to the solid support. This makes the probe 
more available to hybridization and makes the construct more 
stable. If a probe has 5 bases, one may, e.g., use 3 
universal bases as a linker {FIG. 4) . 

15 Determining Sequence from Hybridization Data 

Sequence assembly may be interrupted where ever a 
given overlapping (N-l) mer is duplicated two or more times. 
Then either of the two N-mers differing in the last nucleotide 
may be used in extending the sequence. This branching point 
20 limits unambiguous assembly of sequence. 

Reassembling the sequence of known oligonucleotides 
that hybridize to the target nucleic acid to generate the 
complete sequence of the target nucleic acid may not be 
accomplished in some cases. This is because some information 
25 may be lost if the target nucleic acid is not in fragments of 
appropriate size in relation to the size of oligonucleotide 
that is used for hybridizing. The quantity of information 
lost is proportional to the length of a target being 
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sequenced. However, if sufficiently short targets are used, 

their sequence msy be unambiguously determined. 

The probable frequency of duplicated sequences that 
would interfere with sequence assembly which is distributed 
5 along a certain length of DNA may be calculated. This 

derivation requires the introduction of the definition of a 
parameter having to do with sequence organization: the 
sequence subfragment (SF) . A sequence subf ragment results if 
any part of the sequence of a target nucleic acid starts and 

10 ends with an (N-l)mer that is repeated two or more times 

within the target sequence. Thus, subf ragments are sequences 
generated between two points of branching in the process of 
assembly of the sequences in the method of the invention. The 
sum of all subfragments is longer than the actual target 

15 nucleic acid because of overlapping short ends. Generally, 
subfragments may not be assembled in a linear order without 
additional information since they have shared (N-l)mers at 
their ends and starts. Different numbers of subfragments are 
obtained for each nucleic acid target depending on the number 

20 of its repeated (N-l) mers . The number depends on the value 
of N-l and the length of the target. 

Probability calculations can estimate the 
interrelationship of the two factors. If the ordering of 
positive N-mers is accomplished by using overlapping sequences 
25 of length N-l or at an average distance of A 0/ the N-l of a 
fragment Lf bases long is given by equation one: 
N sf =l+A 0 X K X P ( K, L f ) 



Where K greater than or « 2, and P (K, L f ) represents 
the probability of an N-mer occurring K-times on a fragment L f 
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base long. Also, a computer program that is able to form 
subf ragments from the content of N-mers for any given sequence 
is described below in Example 18. 

The number of subfragments increases with the 
5 increase of lengths of fragments for a given length of probe. 
Obtained subf ragments may not be uniquely ordered among 
themselves. Although not complete, this information is very 
useful for comparative sequence analysis and the recognition 
of functional sequence characteristics. This type of 
10 information may be called partial sequence. Another way of 
obtaining partial sequence is the use of only a subset of 
oligonucleotide probes of a given length. 

There may be relatively good agreement between 
predicted sequence according to theory and a computer 

15 simulation for a random DNA sequence. For instance/ for N-l = 
7, [using an 8-mer or groups of sixteen 10-mers of type 5 1 
(A,T,C,G) B e (A,T,C,G) 3'] a target nucleic acid of 200 bases 
will have an average of three subf ragments . However, because 
of the dispersion around the mean, a library of target nucleic 

20 acid should have inserts of 500 bp so that less than 1 in 2000 
targets have more than three subf ragments . Thus, in an ideal 
case of sequence determination of a long nucleic acid of 
random sequence, a representative library with sufficiently 
short inserts of target nucleic acid may be used. For such 

25 inserts, it is possible to reconstruct the individual target 
by the method of the invention. The entire sequence of a 
large nucleic acid is then obtained by overlapping of the 
defined individual insert sequences. 
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To reduce the need for very short fragments, e.g. 5 0 
bases for 8-mer probes. The information contained in the 
overlapped fragments present in every random DNA fragmentation 
process like cloning, or random PCR is used. It is also 
5 possible to use pools of short physical nucleic acid 

fragments. Using 8-mers or 11-mers like 5' (A, T, C, G) N e (A, 
T, C ,G )3' for sequencing 1 megabase, instead of needing 
20,000 50 bp fragments only 2,100 samples are sufficient. 
This number consists of 700 random 7 kb clones (basic 

10 library) , 1250 pools of 20 clones of 5 00 bp (subf ragments 
ordering library) and 15 0 clones from jumping (or similar) 
library. The developed algorithm (see Example 18) regenerates 
sequence using hybridization data of these described samples. 
ESWl P frE 17 

15 Algorithm 

This example describes an algorithm for generation 
of a long sequence written in a four letter alphabet from 
constituent k- tuple words in a minimal number of separate, 
randomly defined fragments of a starting nucleic acid sequence 

20 where K is the length of an oligonucleotide probe. The 

algorithm is primarily intended for use in the sequencing by 
hybridization (SBH) process. The algorithm is based on 
subf ragments (SF) , informative fragments (IF) and tiie 
possibility of using pools of physical nucleic sequences for 

25 defining informative fragments. 

As described, subfragments may be caused by branch 
points -in the assembly process resulting from the repetition 
of a K-l oligomer sequence in a target nucleic acid. 
Subfragments are sequence fragments found between any two 
30 repetitive words of the length K-l that occur in a sequence. 
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Multiple occurrences of K-l words are the cause of 
interruption of ordering the overlap of K-words in the process 
of sequence generation. Interruption leads to a sequence 
remaining in the form of subfragments. Thus, the unambiguous 
segments between branching points whose order is not uniquely 
determined are called sequence subfragments. 

Informative fragments are defined as fragments of a 
sequence that are determined by the nearest ends of overlapped 
physical sequence fragments. 

A certain number of physical fragments may be pooled 
without losing the possibility of defining informative 
fragments. The total length of randomly pooled fragments 
depends on the length of k- tuples that are used in the 
sequencing process . 

The algorithm consists of two main units. The first 
part is used for generation of subfragments from the set of 
k- tuples contained in a sequence. Subfragments may be 
generated within the coding region of physical nucleic acid 
sequence of certain sizes, or within the informative fragments 
defined within long nucleic acid sequences. Both types of 
fragments are members of the basic library. This algorithm 
does not describe the determination of the content of the 
k-tuples of the informative fragments of the basic library, 
i.e. the step of preparation of informative fragments to be 
used in the sequence generation process. 

The second part of the algorithm determines the 
linear order of obtained subfragments with the purpose of 
regenerating the complete sequence of the nucleic acid 
fragments of the basic library. For this purpose a second, 
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ordering library is used, made of randomly pooled fragments of 
the starting sequence. The algorithm does not include the 
step of combining sequences of basic fragments to regenerate 
an entire, megabase plus sequence. This may be accomplished 
5 using the link-up of fragments of the basic library which is a 
prerequisite for informative fragment generation* 
Alternatively, it may be accomplished after generation of 
sequences of fragments of the basic library by this algorithm, 
using search for their overlap, based on the presence of 
10 common end- sequences . 

The algorithm requires neither knowledge of the 
number of appearances of a given k- tuple in a nucleic acid 
sequence of the basic and ordering libraries, nor does it 
require the information of which k-tuple words are present on 

15 the ends of a fragment. The algorithm operates with the mixed 
content of k-tuples of various length. The concept of the 
algorithm enables operations with the k»tuple sets that 
contain false positive and false negative k- tuples. Only in 
specific cases does the content of the false k-tuples 

20 primarily influence the completeness and correctness of the 
generated sequence . The algorithm may be used for 
optimization of parameters in simulation experiments, as well 
as for sequence generation in the actual SEH experiments e.g. 
generation of the genomic DNA sequence. In optimization of 

25 parameters, the choice of the oligonucleotide probes 

(k-tuples) for practical and convenient fragments and/or the 
choice of the optimal lengths and the number of fragments for 
the defined probes are especially important. 



This part of the algorithm has a central role in the 
30 process of the generation of the sequence from the content of 



WO 98/31836 PCTYUS98/00704 

88 

k-tuples. It is based on the unique ordering of k-tuples by 
means of maximal overlap. The main obstacles in sequence 
generation are specific repeated sequences and false positive 
and/or negative k-tuples. The aim of this part of the 
5 algorithm is to obtain the minimal number of the longest 

possible subf ragments, with correct sequence. This part of the 
algorithm consists of one basic, and several control steps. A 
two-stage process is necessary since certain information can 
be used only after generation of all primary subf ragments . 

10 The main problem of sequence generation is obtaining 

a repeated sequence from word contents that by definition do 
not carry . information on the number of occurrences of the 
particular k-tuples. The concept of the entire algorithm 
depends on the basis on which this problem is solved* In 

15 principle, there are two opposite approaches: 1) repeated 

sequences may be obtained at the beginning, in the process of 
generation of pSFs, or 2) repeated sequences can be obtained 
later, in the process of the final ordering of the 
subf ragments. In the first case, pSFs contain an excess of 

20 sequences and in the second case, they contain a deficit of 
sequences. The first approach requires elimination of the 
excess sequences generated, and the second requires permitting 
multiple use of some of the subfragments in the process of the 
final assembling of the sequence. 

25 The difference in the two approaches in the degree 

of strictness of the rule of unique overlap of k-tuples. The 
less severe rule is: k-tuple X is unambiguously maximally 
overlapped with k-tuple Y if and only if, the rightmost k-l 
end of k-tuple X is present only on the leftmost end of 



WO 98/31836 PCT/US98/00704 

89 

k-tuple Y. This rule allows the generation of repetitive 
sequences and the formation of surplus sequences, 

A stricter rule which is used in the second approach 
has an addition caveat: k-tuple X is unambiguously maximally 
5 overlapped with k-tuple Y if and only if, the rightmost K~l 
end of k-tuple X is present only on the leftmost end of 
k-tuple Y and if the leftmost K-1 end of k- tuple Y is not 
present on the rightmost end of any other k-tuple. The 
algorithm based on the stricter rule is simpler, and is 
10 described herein. 

The process of elongation of a given subfragment is 
stopped when the right k-1 end of the last k-tuple included is 
not present on the left end of any k-tuple or is present on 
two or more k-tuples. If it is present on only one k-tuple the 

15 second part of the rule is tested. If in addition there is a 
k-tuple which differs from the previously included one, the 
assembly of the given subfragment is terminated only on the 
first leftmost position. If this additional k-tuple does not 
exist, the conditions are met for unique k-1 overlap and a 

20 given subfragment is extended to the right by one element . 

Beside the basic rule, a supplementary one is used 
to allow the usage of k-tuples of different lengths. The 
maximal overlap is the length of k-1 of the shorter k-tuple of 
the overlapping pair. Generation of the pSFs is performed 
25 starting from the first k-tuple from the file in which 

k-tuples are displayed randomly and independently from their 
order in a nucleic acid sequence. Thus, the first k-tuple in 
the file is not necessarily on the beginning of the sequence, 
nor on the start of the particular subfragment. The process 
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of subfragment generation is performed by ordering the 
k- tuples by means of unique overlap, which is defined by the 
described rule. Each used k-tuple is erased from the file. 
At the point when there are no further k-tuples unambiguously 
5 overlapping with the last one included, the building of 

subfragment is terminated and the buildup of another pSF is 
started. Since generation of a majority of sub fragments does 
not begin from their actual starts, the formed pSF are added 
to the k-tuple file and are considered as a longer k-tuple. 
10 Another possibility is to form subfragments going in both 

directions from the starting k- tuple. The process ends when 
further overlap, i.e. the extension of any of the 
subfragments, is not possible. 

The pSFs can be divided in three groups: 1) 
15 Subfragments of the maximal length and correct sequence in 
cases of exact k-tuple set; 2) short subfragments, formed 
due to the used of the maximal and unambiguous overlap rule on 
the incomplete set, and/or the set with some false positive 
k~ tuples; and 3) pSFs of an incorrect sequence. The 
20 incompleteness of the set in 2 ) is caused by false negative 
results of a hybridization experiment, as well as by using an 
incorrect set of k-tuples. These are formed due to the false 
positive and false negative k- tuples and can be : a) 
misconnected subfragments; b) subfragments with the wrong 
25 end; and c) false positive k- tuples which appears as false 
minimal subfragments . 

Considering false positive k- tuples, there is the 
possibility for the presence of a k- tuple containing more 
than one wrong base or containing one wrong base somewhere in 
30 the middle, as well as the possibility for a k-tuple with a 
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wrong base on the end. Generation of short, erroneous or 
misconnected subf ragments is caused by the latter k-tuples. 
The k-tuples of the former two kinds represent wrong pSPs with 
length equal to k- tuple length. 

5 In the case of one false negative k- tuple, pSFs are 

generated because of the impossibility of maximal overlapping. 
In the case of the presence of one false positive k-tuple with 
the wrong base on its leftmost or rightmost end, pSFs are 
generated because of the impossibility of unambiguous 
10 overlapping. When both false positive and false negative 

k-tuples with a common k-1 sequence are present in the file, 
pSFs are generated, and one of these pSFs contains the wrong 
k-tuple at the relevant end. 

The process of correcting subfragments with errors 
in sequence and the linking of unambiguously connected pSP is 
performed after subfragment generation and in the process of 
subfragment ordering. The first step which consists of 
cutting the misconnected pSFs and obtaining the final 
subfragments by unambiguous connection of pSFs is described 
below. 

There are two approaches for the formation of 
misconnected subfragments. In the first a mistake occurs when 
an erroneous k-tuple appears on the points of assembly of the 
repeated sequences of lengths k-1. In the second, the repeated 
25 sequences are shorter than k-1. These situations can occur in 
two variants each. In the first variant, one of the repeated 
sequences represents the end of a fragment. In the second 
variant, the repeated sequence occurs at any position within 
the fragment. For the first possibility, the absence of some 
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k- tuples from the file (false negatives) is required to 
generate a mis connect ion . The second possibility requires the 
presence of both false negative and false positive k-tuples in 
the file. Considering the repetitions of k-1 sequence, the 
5 lack of only one k-tuple is sufficient when either end is 

repeated internally. The lack of two is needed for strictly 
internal repetition. The reason is that the end of a sequence 
can be considered inf ormatically as an endless linear array of 
false negative k-tuples. From the "smaller than k-1 case", 
10 only the repeated sequence of the length of k-2, which 

requires two or three specific erroneous k-tuples, will be 
considered. It is very likely that these will be the only 
cases which will be detected in a real experiment, the others 
being much less frequent . 

15 Recognition of the misconnected subfragments is more 

strictly defined when a repeated sequence does not appear at 
the end of the fragment. In this situation, one can detect 
further two subfragments, one of which contains on its 
leftmost, and the other on its rightmost end k-2 sequences 

20 which are also present in the misconnected subf ragment . When 
the repeated sequence is on the end of the fragment, there is 
only one subfragment which contains k-2 sequence causing the 
mistake in subfragment formation on its leftmost or rightmost 
end. 

25 The removal of misconnected subf raments by their 

cutting is performed according to the common rule: If the 
leftmost or rightmost sequence of the length of k™2 of any 
subfragments is present in any other subfragment, the 
subfragment is to be cut into two subfragments, each of them 

30 containing k-2 sequence. This rule does not cover rarer 
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situations of a repeated end when there are more than one 
false negative k- tuple on the point of repeated k-1 sequence. 
Misconnected subf ragments of this kind can be recognized by 
using the information from the overlapped fragments, or 
5 informative fragments of both the basic and ordering 

libraries. In addition, the misconnected subfragment will 
remain when two or more false negative k- tuples occur on both 
positions which contain the identical k-1 sequence. This is a 
very rare situation since it requires at least 4 specific 
10 false k-tuples. An additional rule can be introduced to cut 
these subf ragments on sequences of length k if the given 
sequence can be obtained by combination of sequences shorter 
than k-2 from the end of one subfragment and the start of 
another. 

15 By strict application of the described rule, some 

completeness is lost to ensure the accuracy of the output. 
Some of the subfragments will be cut although they are not 
misconnected since they fit into the pattern of a misconnected 
subfragment. There are several situations of this kind. For 

20 example, a fragment, beside at least two identical k-1 

sequences, contains any k-2 sequence from k-1 or a fragment 
contains k-2 sequence repeated at least twice and at least one 
false negative k-tuple containing given k-2 sequence in the 
middle, etc. 

25 The aim of this part of the algorithm is to reduce 

the number of pSFs to a minimal number of longer subf ragments 
with correct sequence . The generation of unique longer 
subfragments or a complete sequence is possible in two 
situations. The first situation concerns the specific order 

30 of repeated k-1 words. There are cases in which some or all 
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maximally extended pSPs (the first group of pSFs) can be 
uniquely ordered* For example, in fragment 

S~Rl-a-R2-b~Rl-c-R2-E where S and E are the start and end of a 
fragment, a, b , and c are different sequences specific to 
5 respective subfragments and Rl and R2 are two k>l sequences 
that are tandemly repeated, five subf ragments are generated 
(S-Rl, Rl-a-R2, R2-b-Rl, R1-C-R2, and R-E) . They may be 
ordered in two ways; the original sequence above or S-Rl- 
c-R-b-Rl-a-R-E . In contrast, in a fragment with the same 

10 number and types of repeated sequences but ordered 

differently, i.e. S-Rl-a-Rl-b-R-c-R-E, there is no other 
sequence which includes all subf ragments . Examples of this 
type can be recognized only after the process of generation of 
pSFs . They represent the necessity for two steps in the 

15 process of pSF generation. The second situation of generation 
of false short subf ragments on positions of nonrepeated k~l 
sequences when the files contain false negative and /or 
positive k- tuples is more important. 

The solution for both pSF groups consists of two 
20 parts. First, the false positive k- tuples appearing as the 
nonexisting minimal subf ragments are eliminated. All k- tuple 
subfragments of length k which do not have an overlap on 
either end, of the length of longer than k~a on one end and 
longer than k-b on the other end, are eliminated to enable 
25 formation of the maximal number of connections . In our 

experiments, the values for a and b of 2 and 3, respectively, 
appeared to be adequate to eliminate a sufficient number of 
false positive k-tuples. 



The merging of subfragments that can be uniquely 
30 connected is accomplished in the second step. The rule for 



AVO 98/31836 PCT/US98/00704 

95 

connection is: two subfragments may be unambiguously 
connected if, and only if, the overlapping sequence at the 
relevant end or start of two subfragments is not present at 
the start and/or end of any other subf ragment . 

5 The exception is if one subf ragment from the 

considered pair has the identical beginning and end. In that 
case connection is permitted, even if there is another 
subf ragment with the same end present in the file. The main 
problem here is the precise definition of overlapping 

10 sequence. The connection is not permitted if the overlapping 
sequence unique for only one pair of subfragments is shorter 
than k-2, of it is k-2 or longer but an additional subf ragment 
exists with the overlapping sequence of any length longer than 
k>4. Also, both the canonical ends of pSFs and the ends after 

15 omitting one (or few) last bases are considered as the 
overlapping sequences. 

After this step some false positive k-tuples (as 
minimal subfragments) and some subfragments with a wrong end 
may survive. In addition, in very rare occasions where a 

20 certain number of some specific false k-tuples are 

simultaneously present, an erroneous connection may take 
place. These cases will be detected and solved in the 
subf ragment ordering process, and in the additional control 
steps along with the handling of uncut "misconnected" 

25 subfragments. 

The short subfragments that are obtained are of two 
kinds. In the common case, these subfragments may be 
unambiguously connected among themselves because of the 
distribution of repeated k-1 sequences. This may be done 
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after the process of generation of pSFs and is a good example 
of the necessity for two steps in the process of pSF 
generation. In the case of using the file containing false 
positive and/or false negative k-tuples, short pSFs are 
obtained on the sites of nonrepeated k-1 sequences. 
Considering false positive k-tuples, a k- tuple may contain 
more than one wrong base (or containing one wrong base 
somewhere in the middle), as well as k-tuple on the end. 
Generation of short and erroneous (or misconnected) 
subfragments is caused by the latter k-tuples. The k-tuples 
of the former kind represent wrong pSFs with length equal to 
k- tuple length. 

The aim of merging pSF part of the algorithm is the 
reduction of the number of pSFs to the minimal number of 
longer subfragments with the correct sequence. All k-tuple 
subfragments that do not have an overlap on either end, of the 
length of longer than k-a on one, and longer than k-b on the 
other end, are eliminated to enable the maximal number of 
connections. In this way, the majority of false positive 
k-tuples are discarded. The rule for connection is: two 
subfragments can be unambiguously connected if, and only if 
the overlapping sequence of the relevant end or start of two 
subfragments is not present on the start and/or end of any 
other subf ragment , The exception is a subfragment with the 
identical beginning and end. In that case connection is 
permitted, provided that there is another subfragment with the 
same end present in the file. The main problem here is of 
precise definition of overlapping sequence. The presence of 
at least two specific false negative k- tuples on . the points 
of repetition of k-1 or k-2 sequences, as well as combining of 
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the false positive and false negative k- tuples may destroy or 
"mask" some overlapping sequences and can produce an 
unambiguous, but wrong connection of pSFs . To prevent this, 
completeness must be sacrificed on account of exactness: the 
5 connection is not permitted on the end- sequences shorter than 
k-2, and in the presence of an extra overlapping sequence 
longer than k-4. The overlapping sequences are defined from 
the end of the pSFs, or omitting one, or few last bases. 

In the very rare situations, with the presence of a 
10 certain number of some specific false positive and false 

negative k-tuples, some subfragments with the wrong end can 
survive, some false positive k- tuples (as minimal 
subf ragments) can remain, or the erroneous connection can take 
place. These cases are detected and solved in the 
15 subf ragments ordering process, and in the additional control 
steps along with the handling of uncut, misconnected 
subf ragments . 

The process of ordering of subfragments is similar 
to the process of their generation. If one considers 

20 subfragments as longer k-tuples, ordering is performed by 
their unambiguous connection via overlapping ends , The 
informational basis for unambiguous connection is the division 
of subfragments generated in fragments of the basic library 
into groups representing segments of those fragments . The 

25 method is analogous to the biochemical solution of this 

problem based on hybridization with longer oligonucleotides 
with relevant connecting sequence. The connecting sequences 
are generated as subfragments using the k- tuple sets of the 
appropriate segments of basic library fragments. Relevant 

30 segments are defined by the fragments of the ordering library 
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that overlap with the respective fragments of the basic 
library. The shortest segments are informative fragments of 
the ordering library. The longer ones are several neighboring 
informative fragments or total overlapping portions of 
5 fragments corresponding of the ordering and basic libraries . 
In order to decrease the number of separate samples, fragments 
of the ordering library are randomly pooled, and the unique k- 
tuple content is determined. 

By using the large number of fragments in the 

10 ordering library very short segments are generated, thus 
reducing the chance of the multiple appearance of the k-1 
sequences which are the reasons for generation of the 
subf ragments . Furthermore, longer segments, consisting of the 
various regions of the given fragment of the basic library, do 

15 not contain some of the repeated k-l sequences. In every 
segment a connecting sequence (a connecting subf ragment ) is 
generated for a certain pair of the subfragments from the 
given fragment. The process of ordering consists of three 
steps: (1) generation of the k- tuple contents of each 

20 segment; (2) generation of subfragments in each segment; and 
(3) connection of the subfragments of the segments. Primary 
segments are defined as significant intersections and 
differences of k- tuple contents of a given fragment of the 
basic library with the k-tuple contents of the pools of the 

25 ordering library. Secondary (shorter) segments are defined as 
intersections and differences of the k-tuple contents of the 
primary segments . 



30 



There is a problem of accumulating both false 
positive and negative k-tuples in both the differences and 
intersections. The false negative k-tuples from starting 
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sequences accumulate in the intersections (overlapping parts) , 
as well as false positive k-tuples occurring randomly in both 
sequences, but not in the relevant overlapping region. On the 
other hand, the majority of false positives from either of the 
5 starting sequences is not taken up into intersections. This 
is an example of the reduction of experimental errors from 
individual fragments by using information from fragments 
overlapping with them. The false k- tuples accumulate in the 
differences for another reason. The set of false negatives 

10 from the original sequences are enlarged for false positives 
from intersections and the set of false positives for those 
k-tuples which are not included in the intersection by error, 
i.e. are false negative in the intersection. If the starting 
sequences contain 10% false negative data, the primary and 

15 secondary intersections will contain 19% and 28% false 
negative k- tuples, respectively. On the other hand, a 
mathematical expectation of 77 false positives may be 
predicted if the basic fragment and the pools have lengths of 
500 bp and 10,000 bp, respectively. However, there is a 

20 possibility of recovering most of the "lost" k-tuples and of 
eliminating most of the false positive k-tuples. 

First, one has to determine a basic content of the 
k-tuples for a given segment as the intersection of a given 
pair of the k-tuple contents. This is followed by including 

25 all k- tuples of the starting k-tuple contents in the 

intersection, which contain at one end k~l and at the other 
end k-+ sequences which occur at the ends of two k-tuples of 
the basic set. This is done before generation of the 
differences thus preventing the accumulation of false 

30 positives in that process. Following that, the same type of 
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enlargement of k-tuple set is applied to differences with the 
distinction that the borrowing is from the intersections. All 
borrowed k-tuples are eliminated from the intersection files 
as false positives, 

5 The intersection, i.e. a set of common k-tuples, is 

defined for each pair (a basic fragment) X (a pool of ordering 
library) . If the number of k-tuples in the set is significant 
it is enlarged with the false negatives according to the 
described rule. The primary difference set is obtained by 

10 subtracting from a given basic fragment the obtained 

intersection set. The false negative k-tuples are appended to 
the difference set by borrowing from the intersection set 
according to the described rule and, at the same time, removed 
from the intersection set as false positive k-tuples. When 

.15 the basic fragment is longer than the pooled fragments, this 
difference can represent the two separate segments which 
somewhat reduces its utility in further steps* The primary 
segments are all generated intersections and differences of 
pairs (a basic fragment) X (a pool of ordering library) 

20 containing the significant number of k-tuples. K-tuple sets 
of secondary segments are obtained by comparison of k-tuple 
sets of all possible pairs of primary segments. The two 
differences are defined from each pair which produces the 
intersection with the significant number of k-tuples. The 

25 majority of available information from overlapped fragments is 
recovered in this step so that there is little to be gained 
from the third round of forming intersections, and 
differences. 
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(2) Generation of the subfragments of the segments 
is performed identically as described for the fragments of the 
basic library. 

(3) The method of connection of subfragments 

5 consists of sequentially determining the correctly linked 
pairs of subfragments among the subfragments from a given 
basic library fragment which have some overlapped ends. In 
the case of 4 relevant subfragments, two of which contain the 
same beginning and two having the same end, there are 4 

10 different pairs of subfragments that can be connected. In 

general 2 are correct and 2 are wrong. To find correct ones, 
the presence of the connecting sequences of each pair is 
tested in the subfragments generated from all primary and 
secondary segments for a given basic fragment. The length and 

15 the position of the connecting sequence are chosen to avoid 
interference with sequences which occur by chance. They are 
k+2 or longer, and include at least one element 2 beside 
overlapping sequence in both subfragments of a given pair. 
The connection is permitted only if the two connecting 

20 sequences are found and the remaining two do not exist. The 
two linked subfragments replace former subfragments in the 
file and the process is cyclically repeated. 

Repeated sequences are generated in this step. This 
means that some subfragments are included in linked 
25 subfragments more than once. They will be recognized by 

finding the relevant connecting sequence which engages one 
subfragment in connection with two different subfragments. 

The recognition of misconnected subfragments 
generated in the processes of building pSFs and merging pSFs 
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into longer subf ragments is based on testing whether the 
sequences of subfragments from a given basic fragment exist in 
the sequences of subfragments generated in the segments for 
the fragment. The sequences from an incorrectly connected 
5 position will not be found indicating the misconnected 
subfragments. 

Beside the described three steps in ordering of 
subfragments some additional control steps or steps applicable 
to specific sequences will be necessary for the generation of 
10 more complete sequence without mistakes. 

The determination of which subfragment belongs to 
which segment is performed b comparison of contents of 
k-tuples in segments and subfragments. Because of the errors 
in the k-tuple contents (due to the primary error in pools and 
statistical errors due to the frequency of occurrences of 
k-tuples) the exact partitioning of subfragments is 
impossible. Thus, instead of "all or none" partition, the 
chance of coming from the given segment (P(sf,s)) is 
determined for each subfragment. This possibility is the 
function of the lengths of k-tuples, the lengths of 
subfragments, the lengths of fragments of ordering library, 
the size of the pool, and of the percentage of false k-tuples 
in the file: 
P(sf , s)=(Ck-F)/Lsf , 

where Lsf is the length of subfragment, Ck is the number of 
common k-tuples for a given subfragment /segment pair, and F is 
the parameter that includes relations between lengths of 
k-tuples, fragments of basic library, the size of the pool, 
and the error percentage. 



20 
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Subfragments attributed to a particular segment are 
treated as redundant short pSFs and are submitted to a process 
of unambiguous connection. The definition of unambiguous 
connection is slightly different in this case, since it is 
5 based on a probability that subfragments with overlapping 
end(s) belong to the segment considered* Besides, the 
accuracy of unambiguous connection is controlled by following 
the connection of these subfragments in other segments. After 
the connection in different segments, all of the obtained 

10 subfragments are merged together, shorter subfragments 

included within longer ones are eliminated, and the remaining 
ones are submitted to the ordinary connecting process. If the 
sequence is not regenerated completely, the process of 
partition and connection of subfragments is repeated with the 

15 same or less severe criterions of probability of belonging to 
the particular segment, followed by unambiguous connection. 

Using severe criteria for defining unambiguous 
overlap, some information is not used. Instead of a complete 
sequence, several subfragments that define a number of 

20 possibilities for a given fragment are obtained. Using less 
severe criteria an accurate and complete sequence is 
generated. In a certain number of situations, e.g. an 
erroneous connection, it is possible to generate a complete, 
but an incorrect sequence, or to generate "monster" 

25 subfragments with no connection among them . Thus, for each 

fragment of the basic library one obtains: a) several possible 
solutions where one is correct and b) the most probable 
correct solution. Also, in a very small number of cases, due 
to the mistake in the subfragment generation process or due to 

30 the specific ratio of the probabilities of belonging / no 
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unambiguous solution is generated or one, the most probable 
solution. These cases remain as incomplete sequences, or the 
unambiguous solution is obtained by comparing these data with 
other, overlapped fragments of basic library, 

5 The described algorithm was tested on a randomly 

generated, 50 kb sequence, containing 40% GC to simulate the 
GG content of the human genome . In the middle part of this 
sequence were inserted various All, and some other repetitive 
sequences, of a total length of about 4 kb . To simulate an in 
10 vitro SBH experiment, the following operations were performed 
to prepare appropriate data. 

- Positions of sixty 5 kb overlapping "clones" were 
randomly defined, to simulate preparation of a basic library: 

- Positions of one thousand 500 bp "clones" were 
15 randomly determined to simulate making the ordering library* 

These fragments were extracted from the sequence. Random 
pools of 2 0 fragments were made, and k- tuple sets of pools 
were determined and stored on the hard disk. These data are. 
used in the subfragment ordering phase: For the same density 

20 of clones 4 million clones in basic library and 3 million 
clones in ordering library are used for the entire human 
genome. The total number of 7 million clones is several fold 
smaller than the number of clones a few kb long for random 
cloning of almost all of genomic DNA and sequencing by a 

25 gel-based method. 

From the data on the starts and ends of 5 kb 
fragments, 117 "informative fragments" were determined to be 
in the sequence. This was followed by determination of sets 
of overlapping k-tuples of which the single "informative 
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fragment" consist. Only the subset of k- tuples matching a 
predetermined list were used. The list contained 65% 8-mers, 
30% 9-mers, and 5% 10-12 -mers , Processes of generation and 
the ordering of subfragments were performed on these data. 

5 The testing of the algorithm was performed on the 

simulated data in two experiments. The sequence of 50 
informative fragments was regenerated with the 10 0% correct 
data set (over 20,000 bp), and 26 informative fragments (about 
10,000 bp) with 10% false k- tuples (5% positive and 5% 
10 negative ones) . 

In the first experiment, all subfragments were 
correct and in only one out of 50 informative fragments the 
sequence was not completely regenerated but remained in the 
form of 5 subfragments. The analysis of positions of 

15 overlapped fragments of ordering library has shown that they 
lack the information for the unique ordering of the 5 
subfragments. The subfragments may be connected in two ways 
based on overlapping ends, 1-2-3-4-5 and 1-4-3-2-5. The only 
difference is the exchange of positions of subfragments 2 and 

20 4. Since subfragments 2, 3, and 4 are relatively short (total 
of about 100 bp) , the relatively greater chance existed, and 
occurred in this case, that none of the fragments of ordering 
library started or ended in the subfragment 3 region. 

To simulate real sequencing, some false 
25 ("hybridization") data was included as input in a number of 
experiments. In oligomer hybridization experiments, under 
proposed conditions, the only situation producing unreliable 
data is the end mismatch versus full match hybridization. 
Therefore, in simulation only those k- tuples differing in a 
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single element on either end from the real one were considered 
to foe false positives. These "false" sets are made as 
follows. On the original set of a k- tuples of the informative 
fragment, a subset of 5% false positive k- tuples are added. 
5 False positive k- tuples are made by randomly picking a k- tuple 
from the set, copying it and altering a nucleotide on its 
beginning or end. This is followed by subtraction of a subset 
of 5% randomly chosen k-tuples. In this way the statistically 
expected number of the most complicated cases is generated in 
10 which the correct k-tuple is replaced with a k- tuple with the 
wrong base on the end* 

Production of k-tuple sets as described leads to up 
to 10% of false data. This value varies from case to case, 
due to the randomness of choice of k-tuples to be copied, 

15 altered, and erased. Nevertheless, this percentage 3-4 times 
exceeds the amount of unreliable data in real hybridization 
experiments. The introduced error of 10% leads to the two 
fold increase in the number of subfragments both in fragments 
of basic library (basic library informative fragments) and in 

20 segments. About 10% of the final subfragments have a wrong 

base at the end as expected for the k-tuple set which contains 
false positives (see generation of primary subfragments) . 
Neither the cases of misconnection of subfragments nor 
subfragments with the wrong sequence were observed. In 4 

25 informative fragments out of 26 examined in the ordering 

process the complete sequence was not regenerated. In all 4 
cases the sequence was obtained in the form of several longer 
subfragments and several shorter subfragments contained in the 
same segment. This result shows that the algorithmic 
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principles allow working with a large percentage of false 
data . 

The success of the generation of the sequence from 
its k- tuple content may be described in terms of completeness 
5 and accuracy. In the process of generation, two particular 
situations can be defined: 1) Some part of the information is 
missing in the generated sequence, but one knows where the 
ambiguities are and to which type they belong, and 2) the 
regenerated sequence that is obtained does not match the 

10 sequence from which the k- tuple content is generated, but the 
mistake can not be detected. Assuming the algorithm is 
developed to its theoretical limits, as in the use of the 
exact k- tuple sets, only the first situation can take place. 
There the incompleteness results in a certain number of 

15 subfragments that may not be ordered unambiguously and the 
problem of determination of the exact length of monotonous 
sequences, i.e. the number of perfect tandem repeats. 

With false k- tuples, incorrect sequence may be 
generated. The reason for mistakes does not lie in the 

20 shortcomings of the algorithm, but in the fact that a given 

content of k- tuples unambiguously represents the sequence that 
differs from the original one. One may define three classes 
of error, depending on the kind of the false k- tuples present 
in the file. False negative k-tuples (which are not 

25 accompanied with the false positives) produce "deletions" . 
False positive k-tuples are producing "elongations (unequal 
crossing over) " . False positives accompanied with false 
negatives are the reason for generation of "insertions", alone 
or combined with "deletions" . The deletions are produced when 

30 all of the k-tuples (or their majority) between two possible 
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starts of the subfragments are false negatives. Since every 
position in the sequence is defined by k k- tuples, the 
occurrence of the deletions in a common case requires k 
consecutive false negatives. (With 10% of the false negatives 
5 and k=8, this situation takes place after every 108 elements) . 
This situation is extremely infrequent even in mammalian 
genome sequencing using random libraries containing ten genome 
equivalents . 

Elongation of the end of the sequence caused by 

10 false positive k-tuples is the special case of "insertions", 

since the end of the sequence can be considered as the endless 
linear array of false negative k-tuples. One may consider a 
group of false positive k-tuples producing subfragments longer 
than one k-tuple. Situations of this kind may be detected if 

15 subfragments are generated in overlapped fragments, like 
random physical fragments of the ordering library. An 
insertion, or insertion in place of a deletion, can arise as a 
result of specific combinations of false positive and false 
negative k-tuples. In the first case, the number of 

20 consecutive false negatives is smaller than k. Both cases 
require several overlapping false positive k-tuples. The 
insertions and deletions are mostly theoretical possibilities 
without sizable practical repercussions since the requirements 
in the number and specificity of false k-tuples are simply too 

25 high. 

In every other situation of not meeting the 
theoretical requirement of the minimal number an the kind of 
the false positive and/or negatives, mistakes in the k-tuples 
content may produce only the lesser completeness of a 
30 generated sequence. 
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SBH, a sample nucleic acid is sequenced by exposing 
the sample to a support -bound probe of known sequence and a 
labeled probe or probes in solution. Wherever the probes 
ligase is introduced into the mixture of probes and sample, 
5 such that, wherever a support has a bound probe and a labeled 
probe hybridized back to back along the sample, the two probes 
will be chemically linked by the action of the ligase. After 
washing, only chemically linked support -bound and labeled 
probes are detected by the presence of the labeled probe* By 

10 knowing the identity of the support -bound probe at a 

particular location in an array, and the identity of the 
labeled probe, a portion of the sequence of the sample may be 
determined by the presence of a label at a point in an array 
on a Format with a sample of three substrate. And not chances 

15 not working are maximally overlapping sequences of all of the 
ligated probe pairs, the sequence of the sample may be 
reconstructed. Not of the sample to be sequenced may be a 
nucleic acid fragment or oligonucleotide of ten base pairs 
("bp") . The sample is preferably four to one thousand bases 

20 in length. 

The length of the probe is a fragment less than ten 
bases in length, and, preferably, is between four and nine 
bases in length. In this way, arrays of support -bound probes 
may include all oligonucleotides of a given length or may 
25 include only oligonucleotides selected for a particular test. 
Where all oligonucleotides of a given length are used, the 
number of central oligonucleotides may be calculated by 4 N 
where N is the length of the probe. 
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pe-Usina Sequencing Chips 

When ligation is employed in the sequencing process, 
then the ordinary oligonucleotide chip cannot be immediately 
5 reused. The inventor contemplates that this may be overcome 
in various ways * 

One may employ ribonucleotides for the second probe, 
probe P, so that this probe may subsequently be removed by 
RNAse treatment. RNAse treatment may utilize RNAse A an 

10 endoribonuclease that specifically attacks single-stranded RNA 
3 to pyrimidine residues and cleaves the phosphate linkage to 
the adjacent nucleotide. The end products are pyrimidine 3 
phosphates and oligonucleotides with terminal pyrimidine 3 
phosphates . RlSTAse A works in the absence of cof actors and 

15 divalent cations. 

To utilize an RNAse, one would generally incubate 
the chip in any appropriate RNAse -containing buffer , as 
described by Sambrook et al . (1989; incorporated herein by 
reference) . The use of 30-50 ul of RNAse-containing buffer 
20 per 8 x 8 mm or 9 x 9 mm array at 3 7°C for between 10 and 6 0 
minutes is appropriate. One would then wash with 
hybridization buffer. 

Although not widely applicable, one could also use 
the uracil base, as described by Craig et al . (1989) , 
25 incorporated herein by reference, in specific embodiments. 
Destruction of the ligated probe combination, to yield a 
re-usable chip, would be achieved by digestion with the E. 
Coli repair enzyme, uracil -DNA glycosylase which removes 
uracil from DNA. 
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One could also generate a specifically cleavable 
bond between the probes and then cleave the bond after 
detection. For example, this may be achieved by chemical 
ligation as described by Shabarova et al . , (1991) and 
5 Dolinnaya et ai . , (1988), both references being specifically 
incorporated herein by reference. 

Shabarova et al . (1991) describe the condensation of 
oligodeoxyribo nucleotides with cyanogen bromide as a 
condensing agent . In their one step chemical ligation 
10 reaction, the oligonucleotides are heated to 91 a C, slowly 

cooled to 0°C / then 1 ul 10 mM BrCN in acetonitrile is added. 

Dolinnaya et al . (1988) show how to incorporate 
phosphoramidiate and pyrophosphate internucleotide bonds in 
DISTA duplexes. They also use a chemical ligation method for 

15 modification of the sugar phosphate backbone of DNA, with a 
water-soluble carbodiimide (CDI) as a coupling agent. The 
selective cleavage of a phosphoamide bond involves contact 
with 15% CH3COOH for 5 min at 95°C. The selective cleavage of 
a pyrophosphate bond involves contact with a pyridine -water 

20 mixture (9:1) and freshly distilled (CF 3 CO) 2 0. 
E XAMPL E 1 9 

Diagnostics Scoring Known Mutations or Full Gene Re sequencing 

In a simple case, the goal may be to discover 
whether selected, known mutations occur in a DNA segment. Less 
25 than 12 probes may suffice for this purpose, for example, 5 

probes positive for one allele, 5 positive for the other, and 
2 negative for both. Because of the small number of probes to 
be scored per sample, large numbers of samples may be analyzed 
in parallel. For example, with 12 probes in 3 hybridization 
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cycles, 96 different genomic loci or gene segments from 64 
patient may be analyzed on one 6 x 9 in membrane containing 12 
x 24 subarrays each with 64 dots representing the same DNA 
segment from 64 patients. In this example, samples may be 
5 prepared in sixty-four 96-well plates. Each plate may 

represent one patient, and each well may represent one of the 
DNA segments to be analyzed. The samples from 64 plates may be 
spotted in four replicas as four quarters of the same 
membrane . 

10 A set of 12 probes may be selected by single channel 

pipetting or by a single pin transferring device (or by an 
array of individually-controlled pipets or pins) for each of 
the 96 segments, and the selected probes may be arrayed in 
twelve 96-well plates. Probes may be labelled, if they are not 

15 prelabelled, and then probes from four plates may be mixed 
with hybridization buffer and added to the subarrays 
preferentially by a 96 -channel pipeting device* After one 
hybridization cycle it is possible to strip off 

previously-applied probes by incubating the membrane at 37° to 
20 55°C in the preferably undiluted hybridization or washing 
buffer. 

The likelihood that probes positive for one allele 
are positive and probes positive for the other allele are 
negative may be used to determine which of the two alleles is 
25 present. In this redundant scoring scheme, some level (about 
10%) of errors in hybridization of each probe may be 
tolerated . 



An incomplete set of probes may be used for scoring 
most of the alleles, especially if the smaller redundancy is 
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sufficient, e.g. one or two probes which prove the presence or 
absence in a sample of one of the two alleles. For example, 
with a set of four thousand 8-mers there is a 91% chance of 
finding at least one positive probe for one of the two alleles 
5 for a randomly selected locus . The incomplete set of probes 
may be optimized to reflect G+C content and other biases in 
the analyzed samples. 

For full gene sequencing, genes may be amplified in 
an appropriate number of segments. For each segment, a set of 
probes (about one probe per 2-4 bases) may be selected and 
hybridized. These probes may identify whether there is a 
mutation anywhere in the analyzed segments. Segments (i.e., 
subarrays which contain these segments) where one or more 
mutated sites are detected may be hybridized with additional 
probes to find the exact sequence at the mutated sites. If a 
DKTA sample is tested by every second 6-mer, and a mutation is 
localized at the position that is surrounded by positively 
hybridized probes TGCAAA and TATTCC and covered by three 
negative probes: CAAAAC, AAACTA and ACTATT , the mutated 
nucleotides must be A and/or C occurring in the normal 
sequence at that position. They may be changed by a single 
base mutation, or by a one or two nucleotide deletion and/or 
insertion between bases AA, AC or CT . 

One approach is to select a probe that extends the 
25 positively hybridized probe TGCAAA for one nucleotide to the 
right, and which extends the probe TATTCC one nucleotide to 
the left. With these 8 probes (GCAAAA, GCAAAT, GCAAAC, GCAAAG 
and ATATTC , TTATTC , CTATTC , GTATTC) two questionable 
nucleotides are determined. 



15 
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The most likely hypothesis about the mutation may be 
determined. For example, A is found to be mutated to G. There 
are two solutions satisfied by these results. Either 
replacement of A with G is the only change or there is in 
5 addition to that change an insertion of some number of bases 
between newly determined G and the following C. If the result 
with bridging probes is negative these options may then be 
checked first by at least one bridging probe comprising the 
mutated position (AAGCTA) and with an additional 8 probes: 
10 CAAAGA, CAAAGT, CAAAGC, CAAAGG and ACTATT , TCTATT, CCTATT , 
GCTATT , I There are many other ways to select 
mutation- solving probes. 

In the case of diploid, particular comparisons of 
scores for the test samples and homozygotic control may be 
15 performed to identify heterozygotes (see above) . A few 

consecutive probes are expected to have roughly twice smaller 
signals if the segment covered by these probes is mutated on 
one of the two chromosomes . 
EXAMPLE 20 

20 Identification of Genes (Mutations) Responsible 
tor QefoeticPxsoraeyg and Other Traits 



Using universal sets of longer probes (8-mers or 
9-mers) on immobilized arrays of samples, DNA fragments as 
25 long as 5-20 kb may be sequenced without subcloning . 

Furthermore, the speed of sequencing readily may be about 10 
million bp/day/hybridization instrument . This performance 
allows for resequencing a large fraction of human genes or the 
human genome repeatedly from scientifically or medically 
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interesting individuals. To resequence 50% of the human genes, 
about 100 million bp is checked. That may be done in a 
relatively short period of time at an affordable cost. 

This enormous resequencing capability may be used in 
5 several ways to identify mutations and/or genes that encode 

for disorders or any other traits. Basically, mRNAs (which may 
be converted into cDNAs) from particular tissues or genomic 
DNA of patients with particular disorders may be used as 
starting materials. From both sources of DNA, separate genes 

10 or genomic fragments of appropriate length may be prepared 
either by cloning procedures or by in vitro amplification 
procedures (for example by PCR) . If cloning is used, the 
minimal set of clones to be analyzed may be selected from the 
libraries before sequencing. That may be done efficiently by 

15 hybridization of a small number of probes, especially if a 
small number of clones longer than 5 kb is to be sorted. 
Cloning may increase the amount of hybridization data about 
two times, but does not require tens of thousands of PCR 
primers . 

20 In one variant of the procedure, gene or genomic 

fragments may be prepared by restriction cutting with enzymes 
like Hga I which cuts DNA in following way: 
GACGC(N5 1 ) /CTGCG(N10 ' ) . Protruding ends of five bases are 
different for different fragments. One enzyme produces 

25 appropriate fragments for a certain number of genes . By 

cutting cDNA or genomic DNA with several enzymes in separate 
reactions, every gene of interest may be excised 
appropriately. In one approach, the cut DNA is fractionated by 
size. DNA fragments prepared in this way (and optionally 

30 treated with Exonuclease III which individually removes 
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nucleotides from the 3 ' end and increases length and 
specificity of the ends) may be dispensed in the tubes or in 
multiwell plates. From a relatively small set of DNA adapters 
with a common portion and a variable protruding end of 
appropriate length, a pair of adapters may be selected for 
every gene fragment that needs to be amplified. These adapters 
are ligated and then PCR is performed by universal primers. 
From 1000 adapters, a million pairs may be generated, thus a 
million different fragments may be specifically amplified in 
the identical conditions with a universal pair of primers 
complementary to the common end of the adapters . 

If a DNA difference is found to be repeated in 
several patients, and that sequence change is nonsense or can 
change function of the corresponding protein, then the mutated 
gene may be responsible for the disorder. By analyzing a 
significant number of individuals with particular traits, 
functional allelic variations of particular genes could be 
associated by specific traits. 

This approach may be used to eliminate the need for 
very expensive genetic mapping on extensive pedigrees and has 
special value when there is no such genetic data or material. 

ge prjnq Single Nucleotide Polymorphisms in Genetic Mapp ing 

Techniques disclosed in this application are 
appropriate for an efficient identification of genomic 
fragments with single nucleotide polymorphisms (SNUPs) . In 10 
individuals by applying the described sequencing process on a 
large number of genomic fragments of known sequence that may 
be amplified by cloning or by in vitro amplification, a 
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sufficient number of DNA segments with SNUPs may be 
identified. The polymorphic fragments are further used as SNUP 
markers. These markers are either mapped previously (for 
example they represent mapped STSs) or they may be mapped 
5 through the screening procedure described below. 

SNUPs may be scored in every individual from 
relevant families or populations by amplifying markers and 
arraying them in the form of the array of subarrays . Subarrays 
contain the same marker amplified from the analyzed 

10 individuals. For each marker, as in the diagnostics of known 
mutations, a set of 6 or less probes positive for one allele 
and 6 or less probes positive for the other allele may be 
selected and scored. From the significant association of" one 
or a group of the markers with the disorder, chromosomal 

15 position of the responsible gene(s) may be determined. Because 
of the high throughput and low cost, thousands of markers may 
be scored for thousands of individuals. This amount of data 
allows localization of a gene at a resolution level of less 
than one million bp as well as localization of genes involved 

20 in polygenic diseases. Localized genes may be identified by 
sequencing particular regions from relevant normal and 
affected individuals to score a mutation (s) . 

PCR is preferred for amplification of markers from 
genomic DNA. Each of the markers require a specific pair of 
25 primers. The existing markers may be convertible or new 
markers may be defined which may be prepared by cutting 
genomic DNA by Hga I type restriction enzymes, and by ligation 
with a pair of adapters. 
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SNUP markers can be amplified or spotted as pools to 
reduce the number of independent amplification reactions. In 
this case, more probes are scored per one sample. When 4 
markers are pooled and spotted on 12 replica membranes, then 
5 48 probes (12 per marker) may be scored in 4 cycles. 
EXAMPLE 22 

Detection and Verification of Identity of DNA Fragments 

DNA fragments generated by restriction cutting, 
cloning or in vitro amplification (e.g. PCR) frequently may be 

10 identified in a experiment. Identification may be performed 
by verifying the presence of a DNA band of specific size on 
gel electrophoresis. Alternatively, a specific oligonucleotide 
may be prepared and used to verify a DNA sample in question by 
hybridization. The procedure developed here allows for more 

15 efficient identification of a large number of samples without 
preparing a specific oligonucleotide for each fragment. A set 
of positive and negative probes may be selected from the 
universal set for each fragment on the basis of the known 
sequences. Probes that are selected to be positive usually are 

20 able to form one or a few overlapping groups and negative 
probes are spread over the whole insert . 

This technology may be used for identification of 
STSs in the process of their mapping on the YAC clones . Each 
of the STSs may be tested on about 100 YAC clones or pools of 

25 YAC clones. DNAs from these 100 reactions possibly are spotted 
in one subarray. Different STSs may represent consecutive 
subarrays. In several hybridization cycles, a signature may be 
generated for each of the DNA samples, which signature proves 
or disproves existence of the particular STS in the given YAC 

30 clone with necessary confidence. 
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To reduce the number of independent PGR reactions or 
the number of independent samples for spotting, several STSs 
may be amplified simultaneously in a reaction or PCR samples 
may be mixed, respectively. In this case more probes have to 
be scored per one dot . The pooling of STSs is independent of 
pooling YACs and may be used on single YACs or pools of YACs. 
This scheme is especially attractive when several probes 
labelled with different colors are hybridized together* 

In addition to confirmation of the existence of a 
DNA fragment in a sample, the amount of DNA may be estimated 
using intensities of the hybridization of several separate 
probes or one or more pools of probes. By comparing obtained 
intensities with intensities for control samples having a 
known amount of DNA, the quantity of DNA in all spotted 
samples is determined simultaneously. Because only a few 
probes are necessary for identification of a DNA fragment, and 
there are N possible probes that may be used for DNA N bases 
long, this application does not require a large set of probes 
to be sufficient for identification of any DNA segment. From 
one thousand 8-mers, on average about 3 0 full matching probes 
may be selected for a 1000 bp fragment. 
EXAMPLE 23 

Identification of Infectious Disease Organisms and Their* 
VarA.a.ufcs, 

DNA-based tests for the detection of viral, 
bacterial, fungal and other parasitic organisms in patients 
are usually more reliable and less expensive than 
alternatives. The major advantage of DNA tests is to be able 
to identify specific strains and mutants, and eventually be 
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able to apply more effective treatment. Two applications are 
described below. 

The presence of 12 known antibiotic resistance genes 
in bacterial infections may be tested by amplifying these 
5 genes. The amplified products from 128 patients may be spotted 
in two subarrays and 24 subarrays for 12 genes may then be 
repeated four times on a 8 x 12 cm membrane. For each gene, 12 
probes may be selected for positive and negative scoring. 
Hybridizations may be performed in 3 cycles. For these tests , 

10 a much smaller set of probes is most likely to be universal. 
For example, from a set of one thousand 8-mers, on average 3 0 
probes are positive in 1000 bp fragments, and 10 positive 
probes are usually sufficient for a highly reliable 
identification. As described in Example 9, several genes may 

15 be amplified and/or spotted together and the amount of the 

given DNA may be determined. The amount of amplified gene may 
be used as an indicator of the level of infection. 

Another example involves possible sequencing of one 
gene or the whole genome of an HIV virus. Because of rapid 

20 diversification, the virus poses many difficulties for 
selection of an optimal therapy. DNA fragments may be 
amplified from isolated viruses from up to 64 patients and 
resequenced by the described procedure. On the basis of the 
obtained sequence the optimal therapy may be selected. If 

25 there is a mixture of two virus types of which one has the 
basic sequence (similar to the case of heterozygotes) , the 
mutant may be identified by quantitative comparisons of its 
hybridization scores with scores of other samples, especially 
control samples containing the basic virus type only. Scores 

30 twice as small may be obtained for three to four probes that 



WO 98/31836 PCT/US98/00704 

121 

cover the site mutated in one of the two virus types present 
in the sample (see above) . 
EXAMPLE 24 

Forensic anfl Parental identlfigettaflB 

5 Sequence polymorphisms make an individual genomic 

DNA unique. This permits analysis of blood or other body 
fluids or tissues from a crime scene and comparison with 
samples from criminal suspects. A sufficient number of 
polymorphic sites are scored to produce a unique signature of 
10 a sample. SBH may easily score single nucleotide polymorphisms 
to produce such signatures. 

A set of DNA fragments (10-1000) may be amplified 
from samples and suspects. DNAs from samples and suspects 
representing one fragment are spotted in one or several 
15 subarrays and each subarray may be replicated 4 times. In 

three cycles, 12 probes may determine the presence of allele A 
or B in each of the samples, including suspects, for each DNA 
locus. Matching the patterns of samples and suspects may lead 
to discovery of the suspect responsible for the crime. 

20 The same procedure may be applicable to prove or 

disprove the identity of parents of a child. DNA may be 
prepared and polymorphic loci amplified from the child and 
adults; patterns of A or B alleles may be determined by 
hybridization for each. Comparisons of the obtained patterns, 

25 along with positive and negative controls, aide in the 

determination of familial relationships. In this case, only a 
significant portion of the alleles need match with one parent 
for identification. Large numbers of scored loci allow for the 
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avoidance of statistical errors in the procedure or of masking 
effects of de novo mutations. 
EXAMPLE 25 

Assessing Genetic Diversity of Populations or Species 
5 and Biolocrical Diversity of Ecological Niches 



Measuring the frequency of allelic variations on a 
significant number of loci {for example, several genes or 
entire mitochondrial DNA) permits development of different 

10 types of conclusions, such as conclusions regarding the impact 
of the environment on the genotypes, history and evolution of 
a population or its susceptibility to diseases or extinction, 
and others. These assessments may be performed by testing 
specific known alleles or by full resequencing of some loci to 

15 be able to define de novo mutations which may reveal fine 
variations or presence of mutagens in the environment* 

Additionally, biodiversity in the microbial world 
may be surveyed by resequencing evolutionarily conserved DNA 
sequences, such as the genes for ribosomal RNAs or genes for 

20 highly conservative proteins.. DNA may be prepared from the 
environment and particular genes amplified using primers 
corresponding to conservative sequences, DNA fragments may be 
cloned preferentially in a plasmid vector (or diluted to the 
level of one molecule per well in multiwell plates and than 

25 amplified in vitro) . Clones prepared this way may be 

resequenced as described above. Two types of information are 
obtained. First of all, a catalogue of different species may 
be defined as well as the density of the individuals for each 
species. Another segment of information may be used to measure 
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the influence of ecological factors or pollution on the 
ecosystem. It may reveal whether some species are eradicated 
or whether the abundance ratios among species is altered due 
to the pollution. The method also is applicable for 
5 sequencing DNAs from fossils . 
EXAMPLE 2$, 

Detection rvr Quantification of Nucleic Acid Species 

DNA or RNA species may be detected and quantified by 
employing a probe pair including an unlabeled probe fixed to a 

10 substrate and a labeled probe in a solution. The species may 
be detected and quantified by exposure to the unlabeled probe 
in the presence of the labeled probe and ligase. 
Specifically/ the formation of an extended probe by ligation 
of the labeled and unlabeled probe on the sample nucleic acid 

15 backbone is indicative of the presence of the species to be 
detected. Thus, the presence of label at a specific point in 
the array on the substrate after removing unligated labeled 
probe indicates the presence of a sample species while the 
quantity of label indicates the expression level of the 

20 species. 

Alternatively, one or more unlabeled probes may be 
arrayed on a substrate as first members of pairs with one or 
more labeled probes to be introduced in solution. According 
to one method, multiplexing of the label on the array may be 

25 carried out by using dyes which fluoresce at distinguishable 

wavelengths. In this manner, a mixture of cDNAs applied to an 
array with pairs of labeled and unlabeled probes specific for 
species to be identified may be examined for the presence of 
and expression level of cDNA species. According to a 

30 preferred embodiment this approach may be carried out to 
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sequence portions of cDHAs by selecting pairs of unlabeled and 
labeled probes pairs comprising sequences which overlap along 
the sequence of a cDNA to be detected. 

Probes may be selected to detect the presence and 
5 quantity of particular pathogenic organisms genome by 

including in the composition selected probe pairs which appear 
in combination only in target pathogenic genome organisms. 
Thus, while no single probe pair may necessarily be specific 
for the pathogenic organism genome, the combination of pairs 

10 is. Similarly, in detecting or sequencing cDNAs, it might 

occur that a particular probe is not be specific for a cDNA or 
other type of species. Nevertheless, the presence and 
quantity of a particular species may be determined by a result 
wherein a combination of selected probes situated at distinct 

15 array locations is indicative of the presence of a particular 
species . 

An infectious agent with about lOkb or more of DNA 
may be detected using a support-bound detection chip without 
the use of polymerase chain reaction (PCR) or other target 

20 amplification procedures. According to other methods, the 
genomes of infectious agents including bacteria and viruses 
are assayed by amplification of a single target nucleotide 
sequence through PCR and detection of the presence of target 
by hybridization of a labelled probe specific for the target . 

25 sequence. Because such an assay is specific for only a single 
target sequence it therefore is necessary to amplify the gene 
by methods such as PCR to provide sufficient target to provide 
a detectable signal. 
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According to this example, an improved method of 
detecting nucleotide sequences characteristic of infectious 
agents through a Format 3 -type reaction is provided wherein a 
solid phase detection chip is prepared which comprises an 
5 array of multiple different immobilized oligonucleotide probes 
specific for the infectious agent of interest. A single dot 
comprising a mixture of many unlabeled probes complementary to 
the target nucleic acid concentrates the label specific to a 
species at one location thereby improving sensitivity over 
10 diffuse or single probe labeling. Such multiple probes may be 
of overlapping sequences of the target nucleotide sequence but 
may also be non- overlapping sequences as well as non-adjacent. 
Such probes preferably have a length of about 5 to 12 
nucleotides . 

15 A nucleic acid sample exposed to the probe array and 

target sequences present in the sample will hybridize with the 
multiple immobilized probes. A pool of multiple labeled 
probes selected to specifically bind to the target sequences 
adjacent to the immobilized probes is then applied with the 

20 sample to an array of unlabeled oligonucleotide probe 

mixtures. Ligase enzyme is then applied to the chip to ligate 
the adjacent probes on the sample. The detection chip is then 
washed to remove unhybridized and unligated probe and sample 
nucleic acids and the presence of sample nucleic acid may be 

25 determined by the presence or absence of label . This method 
provides reliable sample detection with about a 1000 -fold 
reduction of molarity of the sample agent. 



30 



As a further aspect of the invention, the signal of 
the labelled probes may be amplified by means such as 
providing a common tail to the free probe which itself 
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comprises multiple chromogenic, enzymatic or radioactive 
labels or which is itself susceptible to specific binding by a 
further probe agent which is multiply labelled. In this way, 
a second round of signal amplification may be carried out, 
5 Labeled or unlabeled probes may be used in a second round of 
amplification. In this second round of amplification, a 
lengthy DNA sample with multiple labels may result in an 
increased amplification intensity signal between 10 to 100 
fold which may result in a total signal amplification of 
10 100,000 fold. Through the use of both aspects of this 

example, an intensity signal approximately 100,000 fold may 
give a positive result of probe-DNA ligation without having to 
employ PCR or other amplification procedures* 

According to a further aspect of the invention an 
15 array or super array may be prepared which consists of a 

complete set of probes, for example 40 96 6-mer probes. Arrays 
of this type are universal in a. sense that they can be used 
for detection or partial to complete sequencing of any nucleic 
acid species. Individual spots in an array may contain single 
20 probe species or mixtures of probes, for example N(l-3) B{4-6) 
N{l-3) type of mixtures that are synthesized in the single 
reaction (N represents all four nucleotides , B one specific 
nucleotide and where the associated numbers are a range of 
numbers of bases i.e., 1-3 means 11 from one to three bases".) 
25 These mixtures provide stronger signal for a nucleic acid 

species present at low concentration by collecting signal from 
different parts of the same long nucleic acid species 
molecule. The universal set ■ of probes may be subdivided in 
many subsets which are spotted as unit arrays separated by 
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barriers that prevent spreading of hybridization buffer with 
sample and labeled probe (a) . 

For detection of a nucleic acid' species with a known 
sequence one of more oligonucleotide sequences comprising both 
5 unlabelled fixed and labeled probes in solution may be 

selected. Labeled probes are synthesized or selected from the 
presynthesized complete sets of, for example, 7-mers. The 
labeled probes are added to corresponding unit arrays of fixed 
probes such that a pair of fixed and labeled probes will 
10 adjacently hybridize to the target sequence such that upon 

administration of ligase the probes will be covalently bound. 

If a unit array contains more than one fixed probe 
(as separated spots or within the same spot) that are positive 
in a given nucleic acid species all corresponding labeled 
15 probes may be mixed and added to the same unit array. The 
mixtures of labeled probes are even more important when 
mixtures of nucleic acid species are tested. One example of a 
complex mixture of nucleic acid species are mRNAs in one cell 
or tissue. 

20 According to one embodiment of the invention unit 

arrays of fixed probes allow use of every possible immobilized 
probe with cocktails of a relatively small number of labeled 
probes. More complex cocktails of labeled probes may be used 
if a multiplex labeling scheme is implemented. Preferred 

25 multiplexing methods may use different fluorescent dyes or 
molecular tags that may be separated by mass spectroscopy. 

Alternatively, according to a preferred embodiment 
of the invention, relatively short fixed probes may be 
selected which frequently hybridize to many nucleic acid 
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sequences. Such short probes are used in combination with a 
cocktail of labeled probes which may be prepared such that at 
least one labeled probe corresponds to each of the fixed 
proves. Preferred cocktails are those in which none of the 
5 labeled probes corresponds to more than one fixed probe. 
EXAMPLE 27 

Interrogation of Segments of the HIV Virus with All Possible 
IQ-mers 

In this example of Format III SBH, an array was 

10 generated on nylon membranes (e.g., Gene Screen) of all 

possible bound 5-mers (1024 possible pentamers) . The bound 5- 
mer oligonucleotides were synthesized with 5' tails of 5'- 
TTTTTT-NNN- 3 ' (N = all four bases A, C, G, T, at this step in 
the synthesis equal molar amounts of all four bases are 

15 added) . These oligonucleotides were precisely spotted onto 
the nylon membrane, the spots were allowed to dry, and the 
oligonucleotides were immobilized by treating the dried spots 
with UV light. Oligonucleotide densities of up to 18 
oligonucleotides per square nanometer were obtained using this 

20 method. After the UV treatment, the nylon membranes were 
treated with a detergent containing buffer at 60~80°C. The 
spots of oligonucleotides were gridded in subarrays of 10 by 
10 spots , and each subarray has 64 5-mer spots and 36 control 
spots. 16 subarrays give 1024 5-mers which encompasses all 

25 possible 5-mers. 

The subarrays in the array were partitioned from 
each other by physical barriers, e.g., a hydrophobic strip, 
that allowed each subarray to be hybridized to a sample 
without cross-contamination from adjacent subarrays. In a 
30 preferred embodiment, the hydrophobic strip is made from a 
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solution of silicone (e.g., household silicone glue and seal 
paste) in an appropriate solvent (such solvents are well known 
in the art) . This solution of silicone grease is applied 
between the subarrays to form lines which after the solvent 
evaporates act as hydrophobic strips separating the cells. 

In this Format III example, the free or solution 
(nonbound) 5-mers were synthesized with 3' tails of 5'-NN-3' 
(N - all four bases A, C, G, T) . In this embodiment, the free 
5-mers and the bound 5-mers are combined to produce all 
possible 10-mers for sequencing a known DNA sequence of less 
than 20 kb, 2 0 kb of double stranded DNA is denatured into 4 0 
kb of single- stranded DNA. This 40 kb of ss DNA hybridizes to 
about 4% of all possible 10-mers. This low frequency of 10- 
mer binding and the known target sequence allow the pooling of 
free or solution (nonbound) 5-mers for treatment of each 
subarray, without a loss of sequence information. In a 
preferred embodiment, 16 probes are pooled for each subarray, 
and all possible 5-mers are represented in 64 total pools of 
free 5-mers. Thus, all possible 10-mers may be probed against 
a DNA sample using 1024 subarrays (16 subarrays for each pool 
of free 5-mers) . 

The target DNA in this embodiment represents two-600 
bp segments of the HIV virus. These 600 bp segments are 
represented by pools of 60 overlapping 30-mers (the 30-mers 
overlap each adjacent 3 0 mer by 2 0 nucleotides) . The pools of 
3 0-mers mimic a target DNA that has been treated using 
techniques well known in the art to shear, digest, and/or 
random PCR the target DNA to produce a random pool of very 
small fragments . 
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As described above in the previous Format III 
examples, the free 5 -mars are labeled with radioactive 
isotopes, biotin, fluorescent dyes, etc. The labeled free 5- 
mers are then hybridized along with the bound 5-mers to the 
5 target DNA, and ligated. In a preferred embodiment, 30 0-10 00 
units of ligase are added to the reaction. The hybridization 
conditions were worked out following the teachings of the 
previous examples. Following ligation and removal of the 
target DNA and excess free probe, the array is assayed to 
10 determine the location of labeled probes (using the techniques 
described in the examples above) . 

The known DNA sequence of the target , and the known 
free and bound 5-mers in each subarray, predict which bound 5- 
mers will be ligated to a labeled free 5-mer in each subarray. 
15 The signal from 20 of these predicted dots were lost and 20 

new signals were gained for each change in the target DNA from 
the predicted sequence. The overlapping sequence of the bound 
5-mers in these ten new dots identifies which free, labeled 5- 
mer is bound in each new dot , 

20 Using the described methods, arrays and pools of 

free, labeled 5-mers, the test HIV DNA sequence was probed 
with all possible 10-mers. Using this Format III approach, we 
properly identified the "wild- type" sequence of the segments 
tested, as well as several sequence "mutants" that were 

25 introduced into these segments. 

Sequencing of Repetitive DMA Sequences 



In one embodiment, repetitive DNA sequences in the 
target DNA are sequenced with "spacer oligonucleotides" in a 
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modified Format III approach. Spacer oligonucleotides of 
varying lengths of the repetitive DNA sequence (the repeating 
sequence is identified on a first SBH run) are hybridized to 
the target DNA along with a first known adjoining 
5 oligonucleotide and a second known, or group of possible 

oligonucleotides adjoining the other side of the spacer (known 
from the first SBH run) . When a spacer matching the length of 
the repetitive DNA segment is hybridized to the target, the 
two adjacent oligonucleotides can be ligated to the spacer. 

10 If the first known oligonucleotide is fixed to a substrate, 
and the second known or possible oligonucleotide (s) is 
labeled, a bound ligation product including the labeled second 
known or possible oligonucleotide (s) is formed when a spacer 
of the proper length is hybridized to the target DNA. 

15 EXAMPLE 2 9 

Sequencing Through Branch Points with Format 111 SBH 

In one embodiment, branch points in the target DNA 
are sequenced using a third set of oligonucleotides and a 
modified Format III approach. After a first SBH run, several 

20 branch points may be identified when the sequence is compiled. 
These can be solved by hybridizing oligonucleotide (s) that 
overlap partially with one of the known sequences leading into 
the branch point and then hybridizing to the target an 
additional oligonucleotide that is labeled and corresponds to 

25 one of the sequences that comes out of the branch point . When 
the proper oligonucleotides are hybridized to the target DNA, 
the labeled oligonucleotide can be ligated to the other (s) . 
In a preferred embodiment, a first oligonucleotide that is 
offset by one to several nucleotides from the branch point is 

30 selected (so that it reads into one of the branch sequences) , 



WO 98/31836 PCTYUS98/00704 

132 

a second oligonucleotide reading from the first and into the 
branch point sequence is also selected, and a set of third 
oligonucleotides that correspond to all the possible branch 
sequences with an overlap of the branch point sequence by one 
5 or a few nucleotides (corresponding to the first 

oligonucleotide) is selected. These oligonucleotides are 
hybridized to the target DNA, and only the ' third 
oligonucleotide with the proper branch sequence (that matches 
the branch sequence of the first oligonucleotide) will produce 
10 a ligation product with the first and second oligonucleotides . 

EXAMPLE 3 0 

Multiplexing Probes for A nalyzing a Target Nucleic Acid 

In this Example, sets of probes are labeled with 
different labels so that each probe of a set can be 

15 differentiated from the other probes in the set. Thus, the 
set of probes may be contacted with target nucleic acid in a 
single hybridization reaction without the loss of any probe 
information. In preferred embodiments, the different labels 
are different radioisotopes, or different flourescent labels, 

20 or different EMLs . These sets of probes may be used in either 
Format I, Format II or Format III SBH. 

In Format I SBH, the set of differently lableled 
probes are hybridized to target nucleic acid which is fixed to 
a substrate under conditions that allow differentiation 
25 between perfect matches one base-pair mismatches. Specific 

probes which bind to the target nucleic acid are identified by 
their different labels and perfect matches are determined, at 
least in part, from this binding information. 
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In Format II SBH, the target nucleic acids are 
labeled with different probes and hybridized to arrays of 
probes. Specific target nucleic acids which bind to the 
probes are identified by their different labels and perfect 
5 matches are determined, at least in part, form this binding 
information. 

In Format III SBH, the set of differently labeled 
probes and fixed probes are hybridized to a target nucleic 
acid under conditions that allow perfect matches to be 
10 differentiated from one base-pair mismatches* Labeled probes 
that are adjacent, on the target, to a fixed probe are bound 
to the fixed probe, and these products are detected and 
differentiated by their different labels. 

In a preferred embodiment, the different labels are 
15 EMLs, which can be detected by electron capture mass 

spectrometry (EC-MS) . EMLs may be prepared from a variety of 
backbone molecules, with certain aromatic backbones being 
particularly preferred, e.g., see Xu et al . , J. Chromatog. 
764:95-102 (1997) . The EML is attached to a probe in a 
20 reversible and stable manner, and after the rpobe is 

hybridized to target nucleic acid, the EML is reomved from the 
probe and identified by standard EC-MS (e.g., the EC-MS may be 
done by a gas chromotograph-mass spectrometer) . 

E XAMPLE 3 1 

25 Detection of Low Frequency Targ et Nucle ic Acids 

Format III SBH has sufficient discrimination power 
to identify a sequence that is present in a sample at 1 part 
to 99 parts of a similar sequence that differs by a single 
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nucleotide. Thus, Format III can be used to identify a 
nucleic acid present at a very low concentration in a sample 
of nucleic acids, e.g., a sample derived from blood. 

In one embodiment, the two sequences are for cystic 
5 fibrosis and the sequences differ from each other by a 

deletion of three nucleotides. Probes for the two sequences 
were as follows, probes distinguishing the deletion from wild 
type were fixed to a substrate, and a labeled contiguous probe 
was common to both. Using these targets and probes, the 
10 deletion mutant could be detected with Format III SBH when it 
was present at one part to ninety nine parts of the wild-type. 

Polaroid Apparatus and Method for Analyzing a Target Nucleic 
Acid 

15 An apparatus for analyzing a nucleic acid can be 

constructed with two arrays of nucleic acids, and an optional 
material that prevents the nucleic acids of the two arrays 
from mixing until such mixing is desired. The arrays of the 
apparatus may be supported by a variety of substrates, 

20 including but not limited to, nylon membranes, nitrocellulose 
membranes, or other materials disclosed above* In preferred 
embodiments, one of the substrate is a membrane separated into 
sectors by hydrophobic strips, or a suitable support material 
with wells which may contain a gel or sponge. In this 

25 embodiment, probes are placed on a sector of the membrane, or 
in the well, the gel, or sponge, and a solution (with or 
without target nucleic acids) is added to the membrane or well 
so that the probes are solubilized. The solution with the 
solubilized probes is then allowed to contact the second array 
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of nucleic acids. The nucleic acids may be, but are not 
limited to, oligonucleotide probes, or target nucleic acids, 
and the probes or target nucleic acids may be labeled. The 
nucleic acids may be labeled with any labels conventionally 
5 used in the art, including but not limited to radioisotopes, 
flourscent labels or electrophore mass labels. 

The material which prevents mixing of the nucleic 
acids may be disposed between the two arrays in such a way 
that when the material is removed the nucleic acids of the two 
10 arrays mix together. This material may be in the form of a 
sheet, membrane, or other barrier, and this material may be 
comprised of any material that prevents the mixing of the 
nucleic acids. 

This apparatus may be used in Format I SBH as 
15 follows: a first array of the apparatus has target nucleic 

acids that are fixed to the substrate, and a second array of 
the apparatus has nucleic acid probes that are labeled and can 
be removed to interrogate the target nucleic acid of the first 
array. The two arrays are optionally separated by a sheet of 
20 material that prevents the probes from contacting the target 
nucleic acid, and when this sheet is removed the probes can 
interrogate the target. After appropiate incubation and 
(optionally) washing steps the array of targets may be "read" 
to determine which probes formed perfect matches with the 
25 target. This reading may be automated- or can be done manually 
(e.g., by eye with an autoradiogram) . In Format II SBH, the 
procedure followed would be similar to that described above 
except that the target is labeled and the probes are fixed. 
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Alternatively, the apparatus may be used in Format 
III SBH as follows: two arrays of nucleic acid probes are 
formed, the nucleic acid probes of either or both arrays may 
be labeled, and one of the arrays may be fixed to its 
5 substrate. The two arrays are separated by a sheet of 

material that prevents the probes from mixing. A Format II 
reaction is initiated by adding target nucleic acid and 
removing the sheet allowing the probes to mix with each other 
and the target. Probes which bind to adjacent sites on the 

10 target are bound together (e.g., by base-stacking interactions 
or by covalently joining the backbones) , and the results are 
read to determine which probes bound to the target at adjacent 
sites* When one set of probes is fixed to the substrate, the 
fixed array can be read to determine which probes from the 

15 other array are bound together with the fixed probes. As with 
the above method, this reading may be automated (e.g. , with an 
EL ISA reader) or can be done manually (e.g., by eye with an 
autoradiogram) . 



20 The present invention is not to be limited in scope 

by the exemplified embodiments which are intended as 
illustrations of single aspects of the invention, and 
compositions and methods which are functionally equivalent are 
within the scope of the invention. Indeed, numerous 

25 modifications and variations in the practice of the invention 
are expected to occur to those skilled in the art upon 
consideration of the present preferred embodiments* 
Consequently, the only limitations which should be placed upon 
the scope of the invention are those which appear in the 

30 appended claims. 
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All references cited within the body of the instant 
specification are hereby incorporated by reference in their 
entirety. 
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CLAIMS 

What is claimed is: 



1 1- A method for confirming sequencing results, 

2 comprising the steps of: 

3 obtaining a sequence from a nucleic acid using SBH; 

4 identifying a set of probes that are complementary 

5 and not exactly complementary to the sequence of the nucleic 

6 acid; 

7 hybridizing the probes to the nucleic acid under 

8 conditions that allow the differentiation of perfect matches 

9 from one base mismatches; 

10 confirming that the probes do not form perfect 

11 matches with the nucleic acid. 

1 2. The method of Claim 1, wherein the SBH is 

2 format I SBH. 

1 3. The method of Claim 1, wherein the SBH is 

2 format III SBH. 

1 4. The method of Claim 1, wherein the set of 

2 probes are not exactly complementary to the sequence of the 

3 nucleic acid. 

1 5. A method for confirming sequencing results, 

2 comprising the steps of: 

3 obtaining a sequence from a nucleic acid using SBH; 

4 selecting at least one primer for the nucleic acid; 
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5 sequencing the nucleic acid with the primer using 

6 Sanger - sequenc ing ; 

7 comparing the sequence of the nucleic acid derived 

8 from SBH to the sequence for the nucleic acid derived from 

9 Sanger sequencing . 

1 6 . A method for ordering a plurality of Sf s from a 

2 nucleic acid sequence, comprising the steps of: 

3 obtaining the sequence from the nucleic acid using 

4 SBH; 

5 identifying a plurality of primers from the sequence 

6 of the plurality of Sfs, whereby the primers can initiate a 

7 replication reaction on the nucleic acid that will read 

8 through a branch point; 

9 sequencing the nucleic acid with the primers using 

10 Sanger- sequencing; 

11 comparing the sequence derived by Sanger -sequencing 

12 of the nucleic acid around the branch point to the sequences 

13 of the Sfs, whereby the order of the Sfs is determined. 

1 7. A plurality of probes for analyzing a nucleic 

2 acid wherein the plurality of probes are used to interrogate 

3 the nucleic acid under conditions whereby the plurality of 

4 probes can be differentiated from each other, 

1 8. The probes of Claim 7, wherein the nucleic acid 

2 has a known sequence and probes are labeled with a label . 

1 9. The plurality of probes of Claim 7, wherein the 

2 plurality of probes are labeled with a plurality of different 
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3 labels, whereby the probes can be differentiated from each 

4 other by the different labels attached to the probes. 

1 10. A set of probes for analyzing a nucleic acid, 

2 comprising a plurality of pools of probes wherein each pool is 

3 used to interrogate a nucleic acid, and wherein the plurality 

4 of probes are labeled with a plurality of different labels, 

5 whereby the probes in each pool can be differentiated from 

6 each other by the different labels attached to the probes. 

1 11. The set of probes in Claim 9, wherein the 

2 plurality of different labels are a plurality of different 

3 radioisotopes. 

1 12. The set of probes of Claim 9, wherein the 

■2 plurality of different labels are a plurality of different 

3 flourescent molecules. 

1 13. The set of probes of Claim 9, wherein the 

2 plurality of different labels are a plurality of different 

3 EMLs. 

1 14. The set of probes in Claim 10, wherein the 

2 plurality of different labels are a plurality of different 

3 radioisotopes. 

1 15. The set of probes of Claim 10, wherein the 

2 plurality of different labels are a plurality of different 

3 flourescent molecules. 
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1 16, The set of probes of Claim 10, wherein the 

2 plurality of different labels are a plurality of different 

3 EMLs . 



1 17 . A method for analyzing a nucleic acid, 

2 comprising the steps of: 

3 providing an array of oligonucleotide probes; 

4 introducing a sample nucleic acid to the array; 

5 adding a plurality of labeled probes to the array 

6 under conditions that allow the differentiation of perfect 

7 matches from one base mismatches, ; 

8 adding ligase to the array; 

9 incubating the ligase, labeled probes, sample 

10 nucleic acid and array probes under conditions whereby labeled 

11 probe is ligated to array probes when the labeled probe is 

12 adjacent to the array probe on the sample nucleic acid; and 

13 detecting the labeled probes that have been ligated 

14 to the array. 



1 18. The method of Claim 17, further comprising the 

2 step of removing unligated labeled probe after the incubation 

3 step. 

1 19. The method of Claim 18, wherein the nucleic 

2 acid has a known sequence and plurality of probes are labeled 

3 with a label. 

1 20. The method of Claim 19, wherein the label is 

2 selected from the group consisting of a radioisotope, a 

3 flourescent molecule, and an EML. 
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1 21. The method of Claim 18, wherein the plurality 

2 of probes are labeled with a plurality of different labels, 

3 whereby the probes can be differentiated from each other by 

4 the different labels attached to the probes. 

1 22. The method of Claim 21, wherein the plurality 

2 of different labels are a plurality of different 

3 radioisotopes . 

1 23. The method of Claim 21, wherein the plurality 

2 of different labels are a plurality of different flourescent 

3 molecules. 

1 24. The method of Claim 21, wherein the plurality 

2 of different labels are a plurality of different EMLs , 

1 25. A method for analyzing a plurality of nucleic 

2 acids, comprising the steps of: 

3 obtaining a sample comprising the plurality of 

4 nucleic acids, wherein a target nucleic acid is present at 

5 least in a ratio of one part to ninety nine parts of a nucleic 

6 acid that is homologous to the target and differs by at least 

7 one nucleotide from the target; 

8 selecting a set of probes that will identify the 

9 target nucleic acid; 

10 mixing the sample and the probes under conditions 

11 that allow the differentiation of perfect matches from one 

12 base mismatches; 

13 identifying whether the probes form a perfect match 

14 with a nucleic acid in the sample. 
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1 26. An appartus for analyzing a nucleic acid, 

2 comprising : 

3 a first array of nucleic acids; 

4 a second array of nucleic acids; 

5 a material disposed between the first and second 

6 arrays that prevents the mixing of the nucleic acids in the 

7 first array with the nucleic acids in the second array. 

1 27. The apparatus of Claim 26, wherein the nucleic 

2 acids in the second array are labeled oligonucleotide probes. 

1 28. The apparatus of Claim 27, wherein the nucleic 

2 acids in the first array are a plurality of sample nucleic 

3 acids . 

1 23. A method for analyzing a target nucleic acid, 

2 comprising the steps of: 

3 providing an array of bound probes of known sequence 

4 fixed to a substrate; 

5 providing an array of labeled probes of known 

6 sequence ; 

7 providing a material disposed between the arrays of 

8 bound and. labeled probes that prevents the mixing of the 

9 probes in the bound and labeled probe arrays; 

10 adding the target nucleic acid to the labeled 

1 1 probes ; 

12 removing the material between the bound and labeled 

13 probes so that the labeled probes bound probes and target 

14 nucleic acids are mixed together under conditions that allow 

15 the differentiation of perfect matches from one base 

16 mismatches; 
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17 joining the bound and labeled probes that are 

18 hybridized to adjacent sites in the target nucleic acid; 

19 detecting the labeled probe that has been joined to 

20 the bound probe array* 

1 30. A method for analyzing a target nucleic acid, 

2 comprising the steps of: 

3 providing an array of bound probes of known sequence 

4 fixed to a substrate; 

5 providing an array of labeled probes of known 

6 sequence ; 

7 providing a material disposed between the arrays of 

8 bound and labeled probes that prevents the mixing of the 

9 probes in the bound and labeled probe arrays; 

10 removing the material between the bound and labeled 

11 probes so that the labeled probes and bound probes are mixed 

12 together; 

13 adding the target nucleic acid to the labeled and 

14 bound probes under conditions that allow the differentiation 

15 of perfect matches from one base mismatches; 

16 ligating the bound and labeled probes that are 

17 hybridized to adjacent sites in the target nucleic acid; 

18 detecting the labeled probe that has been ligated to 

19 the bound probe array. 

1 31. A method for analyzing a target nucleic acid, 

2 comprising the steps of: 

3 providing an array of bound probes of known sequence 

4 fixed to a substrate, wherein some of the bound probes are 

5 complementary to a plurality of first portions of the target 

6 nucleic acid; 
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7 providing an array of labeled probes of known 

8 sequence, wherein some of the labeled probes are complementary 

9 to a plurality of second portions of the target nucleic acid 

10 and wherein specific second portions are adjacent to specific 

11 first portions; 

12 providing a material disposed between the arrays of 

13 bound and labeled probes that prevents the mixing of the 

14 probes in the bound and labeled probe arrays; 

15 adding the target nucleic acid to the labeled 

16 probes; 

17 removing the material between the bound and labeled 

18 probes so that the labeled probes bound probes and target 

19 nucleic acids are mixed together under conditions that allow 

20 the differentiation of perfect matches from one base 

21 mismatches; 

22 joining the bound and labeled probes that are bound 

23 at the specific first and second portions in the target 

24 nucleic acid; 

25 detecting the labeled probe joined to the bound 

26 probe array. 

1 32. A method for analyzing a target nucleic acid, 

2 comprising the steps of ; 

3 providing an array of bound probes of known sequence 

4 fixed to a substrate, wherein some of the bound probes are 

5 complementary to a plurality of first portions of the target 

6 nucleic acid; 

7 providing an array of labeled probes of known 

8 sequence, wherein some of the labeled probes are complementary 

9 to a plurality of second portions of the target nucleic acid 
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10 and wherein specific second portions are adjacent to specific 

11 first portions; 

12 providing a material disposed between the arrays of 

13 bound and labeled probes that prevents the mixing of the 

14 probes in the bound and labeled probe arrays; 

15 removing the material between the bound and labeled 

16 probes so that the labeled probes and bound probes are mixed 

17 together; 

18 adding the target nucleic acid to the labeled and 

19 bound probes under conditions that allow the differentiation 

20 of perfect matches from one base mismatches; 

21 joining the bound and labeled probes that are bound 

22 at the specific first and second portions in the target 

23 mac le i c ac id ; 

24 detecting the labeled probe joined to the bound 

25 probe array . 

1 33. A method for analyzing a target nucleic acid, 

2 comprising the steps of: 

3 providing an array of bound target nucleic acids ; 

4 providing an array of labeled probes of known 

5 sequence; 

6 providing a material disposed between the arrays of 

7 bound target and labeled probes that prevents the mixing of 

8 the target nucleic, acid and labeled probes; 

9 removing the material between the bound target and 

10 labeled probes so that the labeled probes and bound target 

11 nucleic acids are mixed together under conditions that allow 

12 the differentiation of perfect matches from one base 

13 mismatches; 
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determining which labeled probes have formed perfect 
matches with the target DNA. 



